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SUMMARY 


The influence of nonequilibrium radiative energy transfer and the 
effect of probe configuration changes on the flow phenomena around a 
Jovian entry body is investigated. The radiating shock-layer flow is 
assumed to be axisymmetric, viscous, laminar and in chemical equili- 
brium. The radiative transfer equations are derived under nonequili- 
brium conditions which include multilevel energy transitions. The 
equilibrium radiative transfer analysis is performed with an existing 
nongray radiation model which accounts for molecular band, atomic line, 
and continuum transitions. The nonequilibrium results are obtained 
with and without ablation injection in the shock layer. The nonequi- 
librium results are found to be greatly influenced by the temperature 
distribution in the shock layer. In the absence of ablative products, 
the convective and radiative heating to the entry body are reduced 
significantly under nonequilibrium conditions. The influence of 
nonequilibrium is found to be greater at higher entry altitudes. With 
coupled ablation and carbon phenolic injection, 16 chemical species are 
used in the ablation layer for radiation absorption. Equilibrium and 
nonequilibrium results are compared under peak heating conditions. For 
the study of the probe shape change effects, the initial body shapes 
considered are 45-degree sphere cone, 35-degree hyperboloid, and 45- 
degree ellipsoid. In all three cases, the results indicate that the 

shock- layer flow field and heat transfer to the body are influenced 
significantly by the probe shape change. The effect of shape change 
on radiative heating of the afterbodies is found to be considerably 
larger for the sphere cone and ellipsoid than for the hyperboloid. 




INTRODUCTION 


During the exploratory missions to the planetary atmospheres, the 
entry spacecraft is subjected to various flow environments and heating 
conditions. The type and intensiveness of this heating mainly depends 
on the atmospheric composition of the planet and the trajectory of the 
entry vehicle. At hypersonic entry conditions, radiation plays a very 
important role in the analysis of the flow phenomena around the planetary 
entry body. The radiative energy transferred to the body from the high- 
temperature shock— layer gas exceeds the convective and aerodynamic 
heating rates. The problem of radiative heating of planetary entry 
bodies has been investigated extensively in the literature (refs. 

1 - 4) . One such situation where the radiative heating constitutes 
the major portion of the heat transferred to the probe is the case of 
Jovian entry heating (refs. 5-8). In order to study the composition 
and structure of Jupiter’s atmosphere, the National Aeronautics and 
Space Administration has scheduled a Jovian mission for a probe space- 
craft in 1985. At times, the entry velocity of the Jupiter probe is 
expected to exceed 40 km/sec. The probe has to survive the intense 
radiative heating during this high-speed entry mission. For this 
purpose, the probe is normally coated with ablative heat shield materi- 
als. As the probe advances through the Jovian atmosphere the heat 
shield ablates and forms a protective layer of cool gases around the 
probe, and these are mainly responsible for absorbing the incoming 
radiation. For the Jovian entry probe, it has been estimated that 
the heat shield weight for thermal protection will be as much as 45 
percent of the total weight of the probe. Since experimental facilities 
cannot adequately simulate the conditions expected during a Jovian 
entry mission, most of the required information must be obtained from 
theoretical studies. This is particularly true for investigating the 
extent of radiative heating to the entry body. 

In order to assess the magnitude of radiative heating to the Jovian 
entry body (and its influence on convective heating and other flow 
phenomena), it is essential to employ meaningful radiative transport 
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models in analyzing the shock-layer flow phenomena. In formulating 
the transfer equations for radiating gases, it is normally assumed 
that the gas is in local thermodynamic equilibrium (LTE) . With this 
assumption, the transfer equation is simplified significantly as the 
populations of the various atomic and molecular states are given by 
the equilibrium (Boltzmann) distribution. For a wide range of condi- 
tions encountered in high-speed gas dynamics, the transitions to or 
from the excited energy states are primarily due to atomic and molecu- 
lar collisional processes (rather than radiative emission and absorption 
processes). Under these conditions, the assumption of LTE is usually 
justified. There are situations, however, where this assumption cannot 
be justified and conditions of nonlocal thermodynamic enuilibrium 
(NLTE) may exist. In a very strict sense, NLTE corresponds to the 
conditions where population densities of various energy levels deviate 
from the equilibrium (Boltzmann) distribution. Quite often it has 
been speculated that use of an LTE radiative transport model may not 
be justified in a shock-heated plasma. This is because, for the com- 
bined conditions of high velocities, high temperatures and low densities, 
the probability of a radiative transition becomes comparable with the 
probability of corresponding collisional transitions (ref. 9). Unless 
the gas is optically thick, the emission of a photon is a process that 
is not balanced by its inverse. Consequently,, the population distri- 
bution among the energy levels departs from that predicted by the 
Boltzmann equation. Only one transition in one atom, molecule or ion 
need be unbalanced in this way to invalidate the LTE assumption. 

Most analyses available on the NLTE radiative heat transfer are 
limited to vibrational ly excited, infrared active diatomic and triatomic 
molecules (refs. 10 - 13). This situation is encountered in some en- 
gineering and upper atmospheric studies. The radiative processes 
associated with the Jovian entry conditions, however, correspond to 
ultraviolet radiation and involve electronic transitions. Furthermore, 
the shock-layer gas consists of molecules and atoms as well as charged 
particles. Hence, the studies available on the NLTE radiative heat 
transfer are not particularly suitable for Jovian entry conditions. 
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One particular study by Horton (ref. 14) estimates the importance of 
nonequilibrium radiation during a hypersonic entry into Jupiter’s 
atmosphere. It is concluded, in this study, that the NLTE effect 
depletes the ablation-layer thickness. This leaves the ablative 
products more transparent, and the absorption of radiation by these 
species is less than the equilibrium value. As a result, the net 
radiative heat flux to the probe surface is more under NLTE conditions. 
Although this conclusion is in general agreement with the kind of 
behavior expected in an ablation-contaminated, nonequilibrium, boundary- 
layer flow, the analysis is not complete. It is important to remem- 
ber that the NLTE phenomenon is among the gas particles themselves 
rather than between the gas and the radiation field. Thus, essentially 
no study is available that treats the problem of NLTE radiative transfer 
under planetary entry conditions in a systematic manner. 

Under NLTE, there are more particles in the higher energy levels 
than normally predicted by the equilibrium theory. This is because 
the particles take considerably longer time to establish a deexcitation 
collision in a nonequilibrium field. The overpopulation of the excited 
energy levels leaves the unexcited levels underpopulated. As a direct 
consequence of this, the absorption pattern of the particles is not the 
same as the equilibrium absorption behavior (since a particle in an 
unexcited state is capable of absorbing more radiation than the one 
in an excited state). Hence, a more detailed analysis of the absorption 
cross section of these particles has to be made under both LTE and NLTE 
conditions. The frequency-dependent absorption coefficient for a non- 
gray gas analysis may be treated either in detail or by a "step model.” 
There are several methods available in the literature for detailed 
computation of the equilibrium absorption coefficient (refs. 15 - 17). 

In a step model, the frequency dependence is broken into a number of 
discrete steps, Sutton (see Zoby et al., ref. 18) developed a 58-step 
model for the hydrogen and helium mixture. For this case, a 30-step 
model developed by Tiwari and Subramanian simplifies the analysis even 
further, and the results compare well with the other 2 models (ref. 19). 
In general, the nonequilibrium absorption coefficient is obtained by 



multiplying the equilibrium values by a nonequilibrium factor. This 
factor is a function of the intensity of the incident radiation, col- 
lisional relaxation time of the particle under investigation, and the 
radiative lifetime of the excited state. It is, therefore, important 
to find an appropriate relation for this nonequilibrium factor for the 
Jovian entry conditions. 

After evaluation of the spectral absorption coefficient, the 
next step is to develop an appropriate expression for spectral and 
total radiative heat fluxes. The expression for the total radiative 
transport involves integration over both the frequency spectrum and 
the physical space. In general, the problem of radiative exchange 
is a complex three-dimensional phenomenon. The inclusion of nonequi- 
librium formulation adds to this complexity with different chemical 
species at different energy levels. Hence, reasonable assumptions are 
required in order to obtain meaningful solutions of the transfer 
equation. In many physically realistic problems, the complexity of 
the three-dimensional radiative transfer can be reduced by introduction 
of the "tangent slab approximation." This approximation treats the 
gas layer as a one-dimensional slab in calculation of the radiative 
transport. Radiation in directions other than normal to either the 
body or shock is neglected. The methods for calculating the divergence 
of the radiative flux and other conservation equations are available 
in references 20 to 22. 

Another problem which arises directly from the radiative heating 
rates is the effect of shape change and mass loss on the flow field 
of a Jovian entry body. For exploratory missions to planets such as 
Mars and Venus, the levels of heating to the entry probe are not severe 
enough to significantly change the mass and shape of the probe. In 
contrast, the large radiative heating rates associated with the Jovian 
entry result in massive ablation of the protective heat shield material 
as the probe advances through the atmosphere (ref, 23). This, in turn, 
results in different probe configurations at different stations along 
the entry trajectory. This change of shape can significantly affect 
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the flow pattern and the heating rate distribution around the probe. 
Sutton et al. (refs. 24 - 26) have presented an inviscid flow-field 
analysis for investigating the effects of the recession of the heat 
shield caused by radiative heating to the Jupiter probe. The initial 
probe configurations considered were hyperboloids and sphere cones 
Important studies available on the effects of mass loss, shape change, 
and real-gas aerodynamics of a Jovian atmospheric probe are discussed 
in a survey article by Walberg et al. (ref. 27). It is clear from the 
review of this article that further studies are needed to investigate 
the effects of the shape changes on heating of the forebody and after- 
body of different Jovian entry probes. 

From the literature survey, it is clear that no work is available 
which considers the influence of nonlocal thermodynamic equilibrium 
radiative heating flow phenomena around a Jovian entry body. A few 
available studies are either inappropriate or provide very little 
information for nonequilibrium analysis under Jovian entry conditions. 
The presence of different ablative species further complicates the 
nonequilibrium analysis. For a correct NLTE analysis, it is essential 
to consider the various excitation and deexcitation probablities of 
all important species in the shock layer. It is important to realize 
that, through an entirely different temperature distribution, the NLTE 
radiative transfer influences significantly the convective heating rate 
as well as other flow properties in the shock layer. Since an accurate 
determination of the total heating rate is essential for the design of 
the heat shield and for assessing the survival of the entry probe, it 
is necessary to investigate the extent of NLTE influence on the entire 
shock- layer flow phenomena in a systematic manner. 

The primary objectives of this study, therefore, are twofold: 
first to investigate the influence of NLTE radiative transfer, and then 
to determine the effects of changes in the probe configurations on the 
entire shock- layer flow phenomena. To accomplish this in a systematic 
manner, the present study is divided into four major areas: (1) sig- 

nificance of radiation models on the flow-field solutions, (2) influence 


of NLTE radiation without ablation injection, (3) importance of NLTE 
radiation with ablation injection, and (4) effect of probe shape change 
on the flow phenomena* 

The basic formulation of the entire problem considered for the 
present study is given in the section titled "Basic Formulation." 
Discussions on the radiative transfer equations, radiation absorption 
models, and radiative flux equations (for both the LTE and NLTE 
conditions) are presented under "Radiative Transport Models." 
Information on collisional relaxation times and radiative lifetimes 
of different species (in their appropriate states) are presented next 
("Radiative Lifetimes and Collisional Processes for the Shock-Layer 
Gases") . Thermodynamic and transport properties for each species 
considered in the shock layer are then giveh ("Thermodynamic and 
Transport Properties"), followed by solution procedures for the NLTE 
radiative flux equations and other shock- layer equations in "Method 
of Solution." The entire results of the study are summarized in 
"Results and Discussion," followed by "Conclusions." 
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BASIC FORMULATION 


r ■ 


Introduction 

The physical model and coordinate system for a Jovian entry body 
are shown in figure 1. In this figure, s is the distance measured 
along the body surface and n is the distance normal to the body 
surface. The flow in the shock layer is considered to be axisymmetric, 
steady, radiating and in chemical equilibrium. Both inviscid as well 
as viscous shock-layer analyses are considered in the present study. 

The basic governing equations (along with the appropriate boundary 
conditions) are presented in this section. 

Inviscid Flow Equations 


( 1 ) 

( 2 ) 

(3) 


(4) 




For the physical model considered, the governing equations for 
inviscid flow are expressed as (refs. 28, 29): 

Continuity : 

(9/3s) [ (r + n cos 0)pu] + (3/3n) (Fcpv) = 0 

s-momentum : 

(u/r) (3u/3s) + v(3u/3n) + (uvk/F) + (1/pF) (3p/3s) = 0 
n-momentum : 

(u/F) (3v/3s) + v(9v/3n) - (u^k/F) + p“l(3p/3n) = 0 

Energy: 

(u/F) (3H/3s) + v(3H/3n) - (u/pF) (3p/3s) 

- (v/p) (3p/3n) + (1/p) (div q^^) = 0 

where 

F = 1 + n< 

div qj^ = (3qj^/3n) + qj.[(K/F) + (cos 6/r)] 
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It should be noted that the above equations are written in nondimen- 
sional form. The quantities used to nondimensionalize these equations 
are defined as: 


s = 


* * 
3 /Rn 


n = n*/RjJ 


* * 

li = u /V 


P = 


* * 
P /P„ 


* * *2 
p = p /(P„ ) 


* *2 
h = h /V 


r = 


* * 

^ /Rx, 


* * * 3 * 

~ ^oo ^ K = K /R 


Viscous Flow Equations 


Basic governing equations for viscous shock layer are obtained 
from the steady-state Navier-Stokes equations by keeping terms up 
to second order in the inverse square root of the Reynolds number, 
e , as (refs . 5 , 30) : 

Continuity: 

C9/9s) [(r + n cos 0)pu] + (9/9n) (T^pv) = 0 

s-momentum: 

p[(u/r) (9u/9s) + v(9u/9n) + (uvk/T)] + r~^(9p/9s) 

= e^j(9/9n) (yij;) + u[(2K7r) + (cos 0/c)]^^) 
n-momentum : 

p[[u/r) (9v/9s) + v(9v/9n) - (u^k/T)] + (ap/9n) = 0 
Energy: 

p[(u/r) (9H/9s) + v(9H/9n)] - v(9p/9n) + pu^(vK/r) 

= £2[(9/9n)$i + ^ 2 ] - { + qj^[(K/r) + (cos 0/r)]) 

where 

V = (9u/9n) - (uk/T) 


(5) 


( 6 ) 


(7) 

( 8 ) 


(9a) 
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= (u/Pt) (3H/9n) - (y/Pr)Xl (9c./8n) - 'E \ 

i=l i=l 

+ (]i/Pr) (Pr - 1) u(9u/9n) - (yKU^/r) 

^2 = [(</P) + (cos 0/c)]<J>i 


In the preceding equations, the total enthalpy H = [ (u^ + v2)/2] + h and 

J. is the mass flux relative to the mass velocity and is given by (refs. 
1 

31, 32): 


= (y/Pr) Le(9c^/9n) 


In equation (10), Le represents the multicomponent Lewis number and, 
in the present analysis, is taken to be a constant for all species 
as 1.10. In addition to the quantities given in equation (5), the 
terms used to nondimensionalize the above set of viscous equations 
are 


* * 

Pr = y C /k 
P 


* * * * 
Le - p C D. ./k 
P ij 


* * 

y = y /y 


Ref 


k = k / (yr, n C ) 
^ Ref p,<» 


* -k 

c = c /c 

P P PL 


Ref 


* * * 

= y (V /c ) 

00 p,o° 


J. = J. 

1 1 


V^Ref 
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In addition to the preceding set of equations for the inviscid and 
viscous shock- layer flow, the species continuity equation and equation 
of state are needed to complete the set of equations. The species 
continuity is given by the expression 


p[(u/r) (9C./9S) + v(9C./9n)] = (eVtc) [0/3n) (r^J.)] 


(9b) 

(9c) 

( 10 ) 


an 


(12a) 
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where 


C12b) 
(13) 

(14a) 
(14b) 

where 

k = 0.65206 - 0.04407 £n(X^„) 

Z = 0.67389 - 0.04637 £n(X,,„) 

HZ 

m = 0.95252 - 0.1447 5,n(X,„) 

HZ 

n = 0.97556 - 0.16149 £n(X,,„) 

HZ 

U = V sin e[l + 0.7476(1 - X,,^)1 
t °o H2 

2 3 

CTU = - 545.37 + 61.608 U - 22459 U + 0.039922 U 

t t t 

4 5 

- 0.00035148 U + 0.0000012361 

t t 

2 3 

CHU = 5.6611 - 0.52661 + 0.020376 - 0.00037861 

4 5 

+ 0.0000034265 - 0.000000012206 

= CTU + 61.2 (1-X^2^ 


? = r + n cos 9 

The equation of state for the gas, in general, can be expressed as 

p = P TCR*/M* C J 
P f 

where ^ represents the specific heat of the gas at the free-stream 
conditions. The expression for the equation of state for a hydrogen/ 
helium mixture is given by Zoby et al . (ref. 33) as: 

T* = C.j,[(p*/1013250)^/(p*/0. 001292)^] 

H* = C^[(p*/1013250)"^/(p*/0. 001292)^] (R*T*/M*) 

H O 


l4 
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= CHU - 0.3167 

and represents the mole fraction of H 2 * 


Free-Stream and Boundary Conditions 

The Jupiter atmosphere mainly consists of hydrogen and helium 
gases. In the past, the nominal composition of the atmosphere was 
assumed to be 85 percent hydrogen and 15 percent helium (ref. 34) . 
For different entry times and altitudes, the free-stream conditions 
are different. The values used in the present study at different 
altitudes are given in table 1. 


Table 1. Free-stream conditions for Jovian entry 


Altitude 
Z, km 

Velocity 

■k 

V , km/sec 

Density 

* O 

p^, kg/m^ 

Pressure 
p^, N/m^ 

Temperature 

* 

T , K 

109 

35.207 

7. 197E-4 

443.02 

163.72 

116 

39.09 

4.360E-4 

245.00 

166.91 

138 

44.04 

2. 127E-4 

132.00 

167.02 


The free-stream enthalpy can be calculated by the relation (ref. 35) : 
H = 1.527 RT 

CO CO 

where R = 8.135 Joules/K-mole is the universal gas constant. The 
number density of hydrogen can be calculated by the ideal gas law 
and the relation can be given as 

N,, = (7.243117 X 10^^) (P /T ) 

2 2 

where X„ is the mole fraction of H„. 

H.-. Z 


(15) 


( 16 ) 
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In order to solve the set of governing equations [eqs. (1) to 
(14)], it is essential to specify appropriate boundary conditions 
at the body surface and shock interface. In all cases, the boundary 
conditions immediately behind the shock wave are calculated by using 
the Rankine-Hugoniot relations. In the viscous analysis, the no-slip 
and no- temperature- jump boundary conditions are used. Consequently, 
the velocities at the surface are 


u(o,n) = u =0 
w 

(17) 

v(o,n) = V =0 
w 

(18) 


The boundary condition given by equation (18) is valid only for 
the case with no mass injection. For this case, the temperature at 
the wall is usually specified as 

T = const. (19) 

w 


For the case of ablation mass injection, the wall temperature is either 

specified or calculated. For the calculated wall temperature conditions, 

the wall temperature is the sublimation temperature of the ablator 

surface. Moreover, the ablation process is assumed to be quasi-steady. 

With these assumptions, the expressions for the coupled mass injection 

rate and the sublimation temperature (for the carbon-phenolic heat 

shield material considered in this study) are given by 

N 


m = 


|(-q -qn )/ [A Cc. h.) -hj|/(p V*) 

' ^c,w ^R,w-^ '^ 1=1 1 i"^w A-* ’ ~ “ 


( 20 ) 


T 

Sub 


5 

E 


A ^ + log p £ A^ . ^ 

• 1 1,1 A ^ ^w " 2,j A 

j = l j=l 


+ 


log p 


w 


5 


E 

j=i 



( 21 ) 
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where p* is the wall pressure in atmospheres and is the ablator 

mass fraction at the wall. The A. . coefficients whose values are 
given in references 5, 30, and 36 to 39 are applicable for a free-stream 
gas composition of 89 percent hydrogen and 11 percent helium by mole 
fraction and for ablator mass fractions of 0.4 to 1.0. If the ablator 
mass fraction at the wall is unity, then these coefficients are valid 
for any free-stream gas mixture. For ablation injection, the elemental 
concentrations at the wall are governed by convection and diffusion 
and are given by 

(8c /8n) - (l/e^) (niSc/u) [Cc.) - (c ) ] =0 

1 W W 1 W 1 - 

where Sc is the Schmidt number (Sc = Pr/Le) and (c^) is the 
elemental mass fraction of the solid ablator material at the surface. 

The heat transferred to the wall due to conduction and mass 
diffusion is referred to as the convective heat flux and is given by 
the expression 

N 

s 

~^c w " £^[(k8T/9n) + (y/Sc) ^ h.(3c^/9n)]^ 


The radiative flux emitted from the wall is given by the relation 


( 0 ) 



e 

w 


*4 


oT 


A value of surface emittance e = 0.8 is used in this study. 

w ^ 


( 22 ) 


(23) 


(24) 
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RADIATIVE TRANSPORT MODELS 


Introduction 

An appropriate expression for the radiative flux is needed 

for the solution of the energy equation presented in the previous 
section. This requires a suitable radiative transport model and a 
meaningful spectral model for variation of the absorption coefficient 
of the gas. In this section appropriate expressions for the spectral 
and total radiative flux are given, and a detailed discussion of 
models for the spectral absorption by the hydrogen-helium gas and 
other important ablative species in the shock layer is presented. 

Radiative Flux Equations 

The equations for radiative transport, in general, are integral 
equations which involve integration over both frequency spectrum and 
physical coordinates. In many physically realistic problems, the 
complexity of the three-dimensional radiative transfer can be reduced 
by introduction of the "tangent slab approximation." This approximation 
treats the gas layer as a one-dimensional slab in calculation of the 
radiative transport. Radiation in directions other than normal to 
either the body or shock is neglected in comparison. Discussions on 
the validity of this approximation for planetary entry conditions are 
given in references 40 to 43. The tangent slab approximation is 
employed in this study. It should be pointed out here that this 
approximation is used only for the radiative transport calculations 
and not for other flow variables. 

LTE radiative flux equations . - For the present study, the 
equations of radiative transport are obtained for a gas confined 
between two infinite, parallel boundaries, the shock wave and the 
body. This is shown in figure 2. For one-dimensional radiation, 
the equations of transfer for a nonscattering medium in local 
thermodynamic equilibrium (refs. 20, 21) are given by 
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C25a) 


yCdl^/dx ) = B (t ) - 

V V V V 

y(dl"'/dT 3 = B (t ) - I 

V V'^ V V 


where 


y = cos 6 



T 

ov 



K (n) dn 


In the above equations, k:^ and B^ represent the frequency-dependent 
linear absorption coefficient and Planck function, respectively. 
Furthermore, it should be noted that and I~ correspond to 

positive and negative values of y, respectively. The boundary 
conditions for equations (25a) and (25b) can be expressed as 

i"'(t ,y) = l'"(0,y), T = 0 

V V V Y) 


I 

V 




I 

p 




T 

V 


T 

ov 


By employing the above conditions, integration of equations (25a) and 
(25b) results in 

T 


I ("T 

V V 



B (t) 


I (t^ ,P) 

V ^ V 


- 1 

X exp[-(x^ - t)/y]y dt 


I (x ,y) exp[-(x - X )/y] 
V ov ov V 


X - 1 

X 

V 


(25b) 

(26a) 

(26b) 

(26c) 


(27a) 

(27b) 

(28a) 

(28b) 
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Figure 2. Physical model and coordinate system for the radiation field. 
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Equations (28a) and (28b) describe the radiation field in terms of the 
temperature field within the medium. The temperature field is expressed 
by the Planck function. The term 1^(0, y) exp (-x^/y) in equation (28a) 
represents the radiant energy that originated at the body surface which 
has been attenuated by the factor exp(-x^yy) as a result of absorption. 
The integral term represents the augmentation of due to gaseous 

emission. A similar explanation goes for equation (28b), with respect 
to the shock surface. 

Referring to figure 2, the spectral radiative flux is expressed 
in terms of intensity of radiation (ref. 21) as 


^Rv 




4tt 


I cos 

V 


0 dfi = 


' 1 


I 


I^Cx^,y)y dy 


(29) 


By noting that and correspond to positive and negative values 

of y, equation (29) can be expressed as 





+ 

y dy 





C30) 


The substitution of values for and from equations (28a) and 

(28b) into equation (30) results in the one-dimensional expression for 
spectral radiative flux (ref. 21) as 


1^(0, y)e 

-[/; 




J E2(t^ - t)dt 


- (t -t ) /y 

1 (t ,-y)e y dy 

V OV 




X ) 

V 



( 31 ) 
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where E (t) 
n 


is the exponential integral function defined by 


E 

n 


(t) 



0 


The expression for the net radiative flux at any location is given by 







dv 


Often, it is desirable to obtain separate relations for total 
radiative flux going towards the body and the bow shock. Upon de- 
noting the radiative flux towards the shock by and towards the 

body by q , equation (33) can be written as 
K 




= qj^Cn) - qj^( 


Jo o 


ir, (t ) dv - / q„ (x ) dv 

‘Rv^ v^ / ^Rv^ v-' 

o 


where 


qR = 2,r 


m 


i;|^(0,y)e"'^v^*' y dy + r B^(t) 
o ^ o 


E 2 (x^ - t) dt dv 


= 2ti I^(T^^,-y)e y dy 

B^(t) E 2 (t - T^) dtj dv 


T 

V 

For diffuse surfaces, 1^(0, y) and independent of 

direction (i.e., independent of y) and may be expressed in terms of 

surface radiosities B. and B„, as 

Iv 2v 


( 32 ) 


(33) 


(34) 


(35a) 


(35b) 
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irl (o,y] = B- . ttI (t ,ji) = B„ , 
V ^ Iv V ov 2v 


Hence, equations (35a) and (35b) are expressed as 




E3(t ) + 
Iv V 


■/' 


B (t) E 2 (t 

V V 


t) dtj 


dv 


qo(^3 = 2 




■/ 


B £ 3(1 - T ) + IT / °^B (t) E 2 (t - T ) dt I dv 

2 v ov V I V V ' 


‘] 


The expressions for surface radiosities appearing in this equation 
(ref. 21 ) are given by 


®lv = "ivl” E3 Ct^^ 3 * w B^Ct) E2(t) dtl 

®2V = Sv ^ 2 E3(t^^) . 

E2(t„^ - t) dtj 


/ 


OV 


IT B (t) 


where and represent the surface reflectance of the body and 

the shock respectively. For nonreflecting surfaces, = p^^ = 0, 
and equations (38a) and (38b) reduce to 


iv 


= TI 


£, B (T ), B^ = IT 
1 V V w 2V 


£o B (T ) 
2 v v^ s 


Sometimes it is convenient to express the radiative flux equations 
in terms of gas emissivities , defined by 

£ = 1 - 2E3(t - t) , £ = 1 - 2E3(t - T ) 

V ■^''V V 25'- 


(36) 


(37) 


(38a) 


(38b) 


(39) 


(40) 
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By noting that 


d£ = -2 E.(t - t) dt, d£ = 2E_(t - x ) dt 

V 2'-v V 2^ V 


(t 
V ^ 


0) = 1 - 2E3(t^) 


+ + 

e, = e 
b w 


V 


T^) = 1 - 2E3CO) 


= 1 - 2(1/2) = 


0 = e (t = T ) 
V V 


e(t = T )=1- 2E3 (t - t ) = e, = e 

V ov ^ ov V b s 


equations (36) and (37) can be written as 


T 

00 £ -I 

,;cn) E3 (t^) * 7. j dej 

*'0 L o -• 


dv 






2B_ E3(t -t)+7t 
2v ov V 


L 


B (e ) de dv 

V V V ' 


3 


If the radiative flux into the slab at the boundaries is neglected, 
then the first right-hand term in equations (41a) and (41b) vanishes 
and the net radiative flux is given by 






B (e ) de dv 

V ^ V 


Depending upon the particular assumptions made in a physical problem, 
use is made of either equations (36) and (37) , (41a) and (41b) , or (42) 
in obtaining the net radiative heat flux. 

For mathematical convenience, exponential integrals often are 
approximated by appropriate exponential functions. There are a few 
standard procedures for doing this, and these are discussed in 


(41a) 


(41b) 


(42) 
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references 20 and 21. It has been demonstrated (e.g. ref. 21) that 
when the exponential integral of third order is approximated by 

2E3(z) = exp(-2z) (43) 

the radiative transport solutions are exact in the optically thin 
limit, and of satisfactory accuracy in the optically thick limit. 

By using equation (43) , approximate expressions for the gas 
emissivities are obtained from equation (40) as 


:^ = 1 - exp[2(t - T^)] 

(44a) 

+ 

= 1 - exp[2(T - t)] 

V V 

(44b) 


Since E (z) = -E ^ (z) , one could obtain the relation for the 
n n-1 

exponential integral of second order by differentiating equation 
(43) as 


E 2 (z) = exp (-2z) (45) 

Use of equations (43) and (45) could be made directly in equations 
(34) , (35a) and (35b) to obtain appropriate relations for the 
radiative heat flux. 

In this study, use of the exponential kernel approximation, as 
given by equations (43) and (45) , is made for the radiative transport 
in the shock layer. Furthermore, the bow shock is considered trans- 
parent, and the free stream is considered cold and transparent. For 
the evaluation of the equilibrium spectral radiative flux, equations 
(36) and (37) are used in the present analysis, 

NLTE radiative flux equations . - In the previous subsection, the 
transfer equation and the resulting radiative flux equations were 
obtained under the condition of local thermodynamic equilibrium. 

For this situation, molecular distribution in the various energy 
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levels is collision dominated, and the emission process depends on 

the local equilibrium temperature alone. The treatment outlined 

in obtaining these equations is of a macroscopic nature. The MLTE 

situation, on the other hand, involves a study of the individual 

molecules of the radiating system since these particles do not obey 

the equilibrium Boltzmann distribution. The derivation of the i^LTE 

transfer equation employs a quantum mechanics treatment. The equation 

of radiative transfer, in general, may be expressed in terms of Einstein 

coefficients A , B and B (refs. 20, 44) as 
nm nm , mn ^ 


dl /dn = 

V 


0 < [1 - 
V '• 


n m 


(B /B )1 
nm mn ' 






^ A 

nm 

B 

mn 

< (N /N ) (B /B ) - 1 
^ m n mn nm 

\ ■ 'J 


where N represents the number density of the nth level, B is 
n ^ mn 

the Einstein coefficient for absorption, B is the Einstein 

^ nm 

coefficient for stimulated emission, and A is the Einstein 

nm 

coefficient for spontaneous emission. The above equation can be 
written in a simplified form as 


dl /dn = p< (J - I ) 


(46) 


(47) 


where J is the NLTE source function and is the reduced 

V) V 

absorption coefficient which includes the effect of induced emission 
(negative absorption) in the medium and is defined as 


K 

V 


< 

V 


[1 - 


(N /N ) 
^ n m 


(B /B )] 
^ nm mn ■' 


(48) 


In this equation is the equilibrium absorption coefficient. 

The NLTE source function is given in terms of the population ratio 
as 
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( 49 ) 


j = A / / b Icn /N ) (B /B ) -ill 
V nm 1 nm r m ^ mn nm'^ ( 

By making use of the relations between the Einstein coefficients, 
the source function can be expressed (refs. 14, 20, 44, 45) as 

= (2hv3/c2)/[(N^/N^) (g^gj - 1] 


where 


A = (2hvVc^) B 
nm nm 


B 

nm 




mn 


In this equation h represents the Planck constant and and 

g^ are corresponding statistical weight factors for the lower and 
upper energy levels (different for different species) , which are 
assumed to be unity in the present analysis. 

From equations (48) and (50) , it is evident that the state 
population ratio N /N has to be known in order to evaluate the 
nonequilibrium absorption coefficient and NLTE source function. 

This is achieved by the method of detailed balancing of various 
transition processes. The three processes involved in the steady- 
state detailed balancing are the induced absorption, induced emission, 
and spontaneous (stimulated) emission. In the induced absorption 
process, a quantum of radiation of appropriate energy and frequency 
is absorbed, and this results in exciting an atom (or molecule) . 

In the induced emission, a quantum of radiation interacts with an 
excited particle to give emission of another quantum of the same 
energy, and thereby the particle reverts to the lower energy state. 

In the spontaneous emission case, an excited particle spontaneously 
emits a quantum of radiation of the appropriate frequency and reverts 
to a lower energy state. These processes may be expressed as 


(50) 

(51a) 

(51b) 
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M + hv ^ 


M* + hv M + 2hv 

M* M + hv (52) 

where M* denotes the excited state of the species M. The statis- 
tical steady-state equation for a particular electronic state is 
given by 

dNi/dt = 0 

This implies that the sums of the radiative and collisional rates into 
and out of the state 1 must be equal. By employing this criterion, the 
state population ratio for any two levels in a multilevel system 
consisting of k levels is expressed (ref. 10) as 


N^/N^ = ^nk ^km,ny^^ ^mk ^kn,m J (53) 

where P is the sum of the radiative term A , and the collision 
nk nk 

term "^he quantity ^ is the probability for all transitions 

from level k to m not involving n such that for k = m, Q, , =1. 

kk,n 

Upon combining equations (50) and (53) , the source function for 
the transition between levels n and m (containing k inter- 
mediate levels) (ref . 46) is given by 


J 

V 



Av 


(b dv 
V 


+ 


* * * * 

£ + 0]/(l + e + n -60) 


(54) 


31 



where 


k 

* k 

^ " C^\m) ^^m,n 

i^n ,m 


(55a) 


(55b) 


* 

£ 


= C /[B 
nm 


mn 


B^(T)] 


(55c) 


6 = c2/2hv3, I* = I /B (T), e* = 0/B (T) (56) 

V V V 

Equations (53) and (54) simplify considerably if the level of 
transitions involved is lower. For example, in a three-level transi- 
tion the relation for the population ratio N /Nt is obtained from 

J. 4. Q 

equation (53) as 


= [Aio + + Cio + C2o(Bi2l^^ + Ci2)]/Qy 

+ Coi + Cq2(A21 + C 21 + B2il^^)]/Q,p 


(57) 


In this equation, represents the total probability for transitions 

from level 2 (i.e., from the third level) and is given by 


Q-p = A21 + B21 + C21 + C20 


(58) 
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It can be shown easily that, under equilibrium conditions, equation (57) 
reduces to the Boltzmann distribution as 


No/Ni = exp(AE/kT) 


( 59 ) 


where AE = Ei - Eq is the difference of energy (in eV) between levels 
1 and 0 . 

Upon combining equations (50) and (57) , a simplified expression 
for the source function (ref. 46) is obtained as 


= <(2hv^/c^)/ [exp(hv/kT) 


1 ] 


} 


^ [1 + ^^^Av^^Av^ B + 62) 


(60) 


where 


== 4tt ^ B dv r61) 

Av *'Av V ^ 

61 = Co 2 (A 21 + C 21 + B 21 I^^)/(Boi Q^) (62a) 

<^2 = [C 20 CB 12 + C 12 ) - Co 2 (A 21 + C 21 + B 21 I^^)]/(Aio (62b) 

In the preceding equations, n = Aio/Cio represents the ratio of the 
collisional deactivation (or relaxation) time and radiative lifetime 
of the first excited state, and quantities 6 j and &2 are the influ- 
ence factors in the NLTE source function arising from the higher level 
energy transitions. 

For a two-level transition, the expression for the source function, 
as given by equation (60), further simplifies (refs. 12 , 47 ) as 
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It is evident from the equation (63) as well as from equation 
(60) that the degree of NLTE depends upon the magnitude of the quantities 
n, and Since values of and 62 are always lower than 

n for all particles involving multilevel energy transitions, the 
extent of nonequilibrium is characterized essentially by the parameter 
n. This implies that the major contribution to the source function 
comes from the transitions (collisional as well as radiative) involving 
the ground state and the first excited state. Consequently, for 
n << 1, the source function becomes the Planck function, and the 
assumption of LTE is justified. In this case, the collisional process 
is sufficiently fast to deexcite the particles to the lower state 
before deexcitation takes place by the emission of radiation. On 
the other hand, the condition of radiative equilibrium is reached 
for n >> 1, and in this case the entire process of excitation and 
deexcitation is radiatively controlled. The NLTE radiation becomes 
important for conditions where n = 0(1). By theoretical considerations, 
Jefferies (ref. 10) has established that the value of <5i^2/P is 
approximately 0.1 for most gases involving multilevel energy transitions. 
Equation (60) is used in evaluating the nonequilibrium source function 
when no ablation mass injection is considered in the study. With 
ablation injection in the shock layer, equation (63) is used for the 
evaluation of to simplify the analysis. 

To find the expression for the NLTE spectral radiative flux, 
the procedure outlined in the previous subsection under the LTE 
conditions [eqs. (29) to (42)] are applicable in general. However, 
in this case the NLTE transfer equation, given by equation (47), 
is integrated between the two parallel boundaries (the body and 
the shock) . The formal solution of equation (47) (ref. 21) is 
given by 
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B, E3CT ) 
Iv v" 


T 

*J Jy(t) E2 (t^ - t) dt 


- - -/'°\ct) EzCt - T^)dt]| 

T ’ 

V 

where is the nonequilibrium optical thickness and is given in 

terms of the nonequilibrium absorption coefficient k* as 

V 

* 

% = / 
o 


Following equation (34), the total radiative flux may be divided 

. + . -- 

into two components as q going towards the shock and q^. going 

R K 

towards the body. Thus equation (64) is rewritten as 


/*°° 

^ Jo •^^iv^ 3 (t^) + 7 t J J^(t)E 2 (x^ - t) dt] dv 


oa 

(n) = 2 f [B £ 3(1 - T ) + Ti r 


ov 


^ J^(t)E 2 (t - T^) dt] dv 

V 


Since NLTE parameter n is the property of the absorbing/emitting 
species, the NLTE source function appears only in the terms which 
are attenuated by absorption of radiation by the gaseous species. 
Equations (64), ( 66 a), and ( 66 b) are used for evaluating the NLTE 
radiative flux. 


Spectral Model for Gaseous Absorption 

Appropriate spectral models for gaseous absorption are needed 
for the solution of the radiative flux equations [eqs. ( 37 ) and ( 66 a)] 
derived under equilibrium and nonequilibrium conditions, respectively. 
The absorption model considered in this study is for a nongray gas 
with molecular band, continuum, and atomic line transitions. In 
general, the spectral absorption coefficient for continuum and line 
transitions nay be expressed as 


(64) 

(65) 

( 66 a) 

( 66 b) 
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( 67 ) 


K = T. k9(v) + T. K^(v) 

V 1 ^ J 

The summations in equation (67) extend over all continuum and 

line transitions, respectively. In the present analysis, only 

+ • 

the transitions of the species H, and H are considered for 
the no ablation case. For the case with carbon phenolic ablation 
injection, the additional chemical species C 2 , C^, C, C'*', C , CO, 

— "f* 

0, 0 ^, e , and 0 are considered for the radiative transport. 

The absorption coefficients for line transitions depend on the 
plasma conditions both through the population of the absorption levels 
and the shape of the spectral lines. For heavy atomic species at high 
temperatures, the dominant mechanism for the line broadening is the 
Stark broadening by electron impacts. Following Armstrong et al. 

(ref. 48), the lines can be treated as having the Lorentz shape, for 
which the shape factor is given by 

bk(vj = CYk/")/[('' - + (Yk3 ] (683 

where is the frequency of the kth line center and is the 

Stark half-width of the line. In calculating the absorption co- 
efficients due to atomic line transitions, a line grouping technique 
is used. In this technique, line transitions near a specified 
frequency value are grouped together, and the spectral absorption 
is given as that from the line group. However, each line within 
the group is treated individually. 

The continuum contribution depends mainly on the plasma state 
through the population of absorbing levels. The spectral absorption 
coefficient due to continuum transitions is given by 

K^(v) = Z N. . a^. (v) (69) 

1 ^ ij ij 
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c 

where N. . is the number density of the absorbing level and a. . is 

ij ij 

its cross section. The number density of the particular particle is 
obtained from thermodynamic state calculations. 

By employing the detailed information oh line and continuum 
absorption, Nicolet developed a fairly sophisticated radiative trans- 
port model for applications to planetary entry environment (refs. 17, 
49, 50) . For calculation of the equilibrium radiative transport 
properties in the ablation contaminated boundary layers, the method 
given in reference 49 is used. An approximate model for the fre- 
quency-dependent absorption coefficient is also developed by Sutton 
(see Zoby et al., ref. 18). In this model, the frequency dependence 
of the absorption coefficient is represented by a step-function 
with 58 steps of fixed (but not necessarily equal) widths. In this 
model, the absorption of helium species is neglected. In step- 
function models, the total absorption coefficient of the jth step is 
a summation of the average absorption coefficient for the ith 
transition in the jth step, given by 


X . (w) - I K . . (w) 

J ij 


K..(co) = ( 1 /Av ) J J J K. (oj) dv 

J YI . ^ 

J 

K. = f(T,N.,v) 

Once again k is the equilibrium absorption coefficient, v the 
frequency in eV, T is the temperature in degrees, and is 

the number density in cm” ^ . In this model, the absorption coef- 
ficients for the free-free and bound-free transitions of atomic 
hydrogen are expressed by 

H "35 1/2 . 

= (2.61 X 10 N ^)/(T v3) 

ff e 


(70a) 

(70b) 

(70c) 

(71) 
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= [a-99 " N^)/v3] 


{ 


expU-157780/T) [1 * (1/ 




( 72 ) 


and 


for 1 ^ j< 4 


H "^0 

= [6.31 X 10 (T/v^) N^] exp(A) [exp(B) - 1] 


(73) 


for 4 < N- < N- 

I — £max 


where 


A = (-157880/T) [1 - (6/13.6)] 


B = (157780/T) [(1/25) - (6/13.6)] 


The reduction in the ionization potential 6 is calculated by 


5 = (1.79 X 10" 5 N^2/7y 


(74) 


For bound-bound transition of hydrogen molecules and atoms 


^bb = S b^Cv) 


(75) 


where the line strength S is given by 


S = (1.10 X iol6 f n| N^) exp|(-157780/T) [1 - (l/n|)]| 


(76) 
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The line-shape factor b (v) is given by equation (68). Using 

1C 

the expressions given above, the absorptions due to continuum and line 
transitions over each step are calculated individually. The total 
absorption over each individual step is a strong function of temper- 
ature, and this model is valid for a wide range of temperatures. 
Further information on the 58- step model is available in reference 
18. This model is used in evaluating the equilibrium absorption 
coefficient for the no ablation injection case. 

At high temperatures, the frequency dependence of the absorption 
coefficient is more orderly because of the relative importance of 
continuum transitions over line transitions. Under such conditions, 
it is possible to represent the spectral absorption of the gas by a 
relatively fewer number of steps. A spectral model consisting of 
30 steps is introduced in this study to represent the absorption by 
the hydrogen species in the spectral range of 0 to 20 eV. The ab- 
sorption by the helium species is also neglected in this study. 

The procedure for developing this model is to calculate the spectral 
absorption coefficient first by employing Nicolet’s detailed model. 

The 30- step model is illustrated in figure 3. Further details on this 
30-step model are given in reference 19. Some results obtained using 
the various models are presented under ''Results and Discussion" for 
comparative purposes. 

The methods outlined thus far to evaluate the equilibrium absorp- 
tion coefficients are applicable for hydrogen atoms and molecules. 

This is quite sufficient if no ablation mass injection is considered 
in the analysis. However, with ablation products in the shock layer, 
appropriate relations are needed to evaluate the continuum and line 
contributions to the absorption coefficient by these species. 
Information for obtaining the absorption cross section for the 
electronic band systems belonging to diatomic molecules is available 
in references 51 and 52. According to these references, the equili- 
brium absorption coefficient for a transition from level 0 to level 1 
(higher level) is given by 






( 77 ) 


K (0,1)= o (0,1) n 
V V o 

where is the number density of the absorbing species in level 0 

and a^(0,l) is the absorption cross section. This equilibrium 
absorption coefficient has to be corrected for the NLTE transition. 

The NLTE absorption coefficient is given by reference 47. 

k^(NLTE) = K^«^l - Bio + Aio Cio)/(Coi + Boi)]j- 

where is the frequency-dependent equilibrium absorption coefficient. 

By noting that 

Aio = (2hv3/c2) Bio 

Bio “ (So/gl) Boi 

Cio (Bio/Boi) Coi exp(hv/kT) 

= 1/Cio, = 1/Aio 

equation (78) becomes 

k*CNLTE) = >c^|l - |[1 + n (n/S) / dv ] 

>[(n/6) + exp(-hv/kT)] ^ | 

This equation is used in evaluating the nonequilibrium absorption 
coefficients . 


(78) 

(79a) 

(79b) 

(79c) 

(79d) 


(80) 


It is seen from equations (60), (63), and (80) that, for the 
evaluation of the nonequilibrium absorption coefficients and source 
function (and hence the nonequilibrium spectral radiative heat flux). 
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one must have information on- the collisional relaxation time and the 
radiative lifetime of the excited states. The procedure for obtaining 
these is discussed in the next section. 
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RADIATIVE LIFETIMES AND COLLISIONAL PROCESSES FOR THE SHOCK- LAYER GASES 


Introduction 

For a gaseous medium in local thermodynamic equilibrium, infor- 
mation on the collisions between the various particles and their 
radiative lifetimes in the excited states is of little importance. 

To analyze the nonequilibrium phenomena, however, a quantitative 
study of collisions between various particles is necessary. This 
study need be extended only to those species whose collisional 
transition rates (under a given set of physical conditions) are 
comparable to their radiative transition rates. Monequilibrium 
phenomena are important only under these conditions. 

For Jovian entry conditions and for the case with no ablation 
from the probe surface, the NLTE effects are considered only for 
the hydrogen species in the shock layer. However, for the case of 
carbon-phenolic ablation injection, the molecules play a very 
important role in the radiation blockage. For this case, therefore, 
contributions of the molecules are also included in the NLTE 
analysis. A discussion of the radiative lifetimes and collisional 
relaxation times of various species considered for the NLTE analysis 
(in the shock- layer gas) is presented in this section. 

Radiative Lifetimes of Excited States 

The radiative lifetime n of an excited state is the inverse 

r 

of the Einstein coefficient for spontaneous emission A^q, and this 
is given (refs. 53, 54) by 

Dp = 1/Aio = (c^/8irhv3)/Boi (81) 

where is the Einstein coefficient for absorption. By using 

the statistical relation for Bqj , equation (81) is written as 
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I. 



(82) 


n 


-1 

r 


(8ttv^/ c^) 


/ 


K dv 
V 


For the present study it is more appropriate to use the radiative 
lifetime in a different form. For a fully allowed electronic transition 
in the visible or near-ultraviolet region of the spectrum, the radiative 
lifetime is expressed in terms of the half-width of the spectral line 
(ref. 55) as 


= h/(27Tb) 


(83) 


Upon inserting the numerical value for the Planck constant h (and 
converting from ergs to cm ^) , the value of is found to be 

= (5.3 X 10-^2^/5 


For a Stark-broadened hydrogen line of Lorentz shape, the value of 
the radiative lifetime is found to be 0.52 x io“^ sec (ref 55). 

This value is used in the present study. It should be noted that 
the radiative lifetime is a function of frequency and, therefore, 
varies according to the various levels of energy transitions. 
Considering all radiative transition probabilities between different 
energy levels, it is found that the above-mentioned value of 
corresponds to the shortest time between the ground level (with prin- 
cipal quantum number 2). Since the maximum nonequilibrium effect 
would correspond to this value of its use is justified in 

the present study. 


There are eight known C 2 -band systems in the 0.1- to 7.0-eV 
(0.2- to 1.2-u) spectral region. The Swan band system, whose 
electronic transition is represented by "th® 

notation of ref. 55), is the strongest radiating system of the C., 


molecules. The Freymark band (E’E - A'lr ) and the ^4ulliken band 

+ + ^ 

(D'E^ - x'E ) are the next important radiating systems of the 
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molecules. Hence, the radiative lifetime of these band systems has 
to be determined first in order to evaluate their respective NLTE 
cross sections. For a molecule with an electronic transition, the 
radiative lifetime is related to the wavelength- dependent electronic 
F-number, the electronic transition moment, and the Einstein coefficient- 
for spontaneous emission The measured value for the electronic 

transition moment is available for the band systems (ref. 56) . 

The electronic F-number is given in terms of the square of the tran- 
sition moment [r /ea^ p as 
' e O' 

^ 2 

F T C^) = [(STT^m c/3he^X) z|r /ea | ] /A (35) 

e 1 e e u 

where 

A = (2 - (2S" + 1) 

In the above discussions and relations, the superscript ’ denotes 
the upper state and ” the lower state, and S" represents the 
spin quantum number of the lower state. The quantity ~ ^ For 

A = 0 and ^ For A 7 ^ 0, and A is the resultant angular 

momentum of electrons. The electronic F-number is related to 
the band oscillator strength by 

^v'v" ^el %'v” (36) 

where q , ,, is the Franck-Condon factor. Now, the radiative lifetime 
Ly t V 

of the v' state can be expressed as 

(np = A^,^„ = F^,^„ [(6.67 X 10'5 )/a 2] (g"/g ' ) (87) 

Here, g' and g" are the degeneracies of the upper and lower levels 
respectively and g"/g' = !• The value of F measured for the 

C^-Swan (0,0) band is 6.5 x 10“^ (ref. 57). Upon substituting this 
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value of using equation (87), the radiative lifetime of 

the C^-Swan 3 tt (v’ = 0) state is found to be 8.0 x 10“^ sec. No 

measured value of the band oscillator strength Pyfytt available 
for the Freymark and Mul liken band system. However the values of 
the electronic transition moment are available for these bands; and, 
since F is directly proportional to the electronic transition 

moment, the radiative lifetime of the state of the Freymark 

(0,1) band is 1.5 x 10"^ sec and the Mulliken (0,0) 

band is 6.95 x 10“^ sec. These values are used in the majority of 
cases investigated in the present study. However, another set of 
values for the radiative lifetime of different bands is suggested 
in the literature (ref. 56). These are ri^(Swan) = 1.25 x 10~^ sec, 
ri^(Freymark) = 4.67 x lO ® sec, and (Mulliken) = 8.77 x 10”^ sec. 
These values are significantly different from those mentioned 
earlier. Because of this discrepancy, it is essential to inves- 
tigate the influence of a different set of radiative lifetimes on 
the NLTE results. The following format, therefore, will be adapted 
for calculating the NLTE results in the presence of the ablative 
products in the shock layer: 


NLTE (I): based on n^(Swan) = 8 x lo ^ sec 

NLTE (II): based on n^(Swan) = 1.25 x 10~^ sec 


NLTE (III): based on the combined contributions of ri^(Swan) 
= 8 X 10“^, (Freymark) = 1.5 x 10“^, and 
n^(Mulliken) = 6.95 x 10~® sec 

NLTE (IV): based on the combined contributions of n^(Swan) 

= 1.25 X 10”^, (Freymark) = 4.67 x 10“®, and 
ri^(Mulliken) = 8.77 x 10“^ sec 
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Collisional Relaxation Time 


The collisional relaxation time is a strong function of temperature 
and pressure. At very high temperatures (associated with the ionization 
mechanism of hydrogen atoms), it is reasonable to assume that most of 
the NLTE transfer of radiation takes place before the ionization of 
the hydrogen atom. Hence, for the case with no ablation from the probe 
surface, it is quite sufficient to consider only the collisional process 
between ^^id H and H^. Collisions between two atoms or molecules 
may be of the first or second kind. In collisions of the first kind, 
the kinetic energy (KE) of translation goes into excitation energy 
according to the process 

A+B + KE-^A+B* (88) 

where A and B are two different (or same) species in the ground 

'k 

state and B is the species B in an excited state. Collisions of 
the second kind are more important for the NLTE analysis at high 
temperatures. This process is described by 

A+B*-^A+B + KE (89) 

Here, an atom or molecule gives up excitation energy by colliding 
with another partner. The time taken by a particle to reach the 
ground state from an excited state is the collisional relaxation 
time (ref. 58) . 

The collisional relaxation time for collisions between neutral 
particles (such as atoms and molecules) is given, in general, by 
the relation 

- 1 
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where is the frequency of collisions (sec~^), n is the number 

density of the colliding particles (cm~^), is the collisional 

cross section of the colliding particles (cm^) , and v is the most 
probable velocity of the particles (cm-sec"^). By making use of the 
gas kinetic relations, the relaxation time (in sec) is expressed in 
terms of temperature as 


^ 1/2 
= 4/2 (7im/8RT) (1/n) 


(91) 


where m is the mass of the colliding particles. Information on the 
collisional cross section for different colliding particles is available 
in references 9 and 59. 

For an extensive study of the NLTE process, a detailed study of 
the various collisional processes between the particles (present in 
the radiation field) is necessary. All collisions which are effective 
in deactivating the excited particles must be considered in the analysis. 
Relaxation times for the case of collisions between neutral particles 
of the same kind (i.e. for collisions) can be 

calculated using equation (91) with the appropriate values for colli- 
sional cross section and mass. On the other hand, the collisions 
between unlike particles with different masses and cross sections 
may be equally effective in removing the electronic excitation energy 
from these particles. In such cases, the effective collisional 
frequency is calculated by the method explained below. The colli- 
sions between H and H 2 are taken as an illustrative example; the 
procedure is the same for other combinations of molecules and atoms. 

The number of deactivation collisions made per second between H and 
is given by 


Z(H - H^) 


2nn [2kT(n.„ 




1/2 


(92) 
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where n„ is the number density of and m and m are the 
1^2 ^ 

masses of a hydrogen molecule and atom respectively. The quantity 

0 , is the effective collisional cross section and is given by the 
E 

expression 

fig = [Cfi„ + fij^)/4] + [(fi fi^)V2/2^] C93) 

where q,, and are the cross sections of and H respectively. 

Hence, the relaxation time (in sec) for collisions is the 

reciprocal of the collisional frequency Z (II 2 -H) . 

The combined relaxation time for self-collisions and collisions 
of different kind is given by reference 60. 

l/n^(C0M) = X/[n^(H-H)] + (1 -X)/[ti^(H2-H)] (94) 

where X represents the mole fraction of hydrogen atoms. In the 
present study, collisions between H-H and are considered in 

evaluating the combined relaxation time with no ablation injection 
in the shock-layer. Collisions between ^ 2~^2 neglected because 

the number density of is small compared to the number density of 
atomic hydrogen. Various collisional relaxation times obtained for 
different collisional processes (between the shock-layer species 
without ablative products) are shown in figure 4 as a function of 
temperature. With the ablation injection, the important species 
influencing the collisional deactivation process of the molecules 
are H, H^, C^, C, and e . The relaxation times obtained by 
using equations (91) to (94) for the H-H and combined colli- 

sions are shown in figure 5 as a function of temperature. The 

^ 2 ”^’ ^2"^^^ collisions represent the combined colli- 

sional process. The radiative lifetimes of the Swan (0,0), Freymark 
(0,1) and Mulliken (0,0) bands, corresponding to the case of (III), 
are also shown in this figure. 
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Figure 5. Variation of collisional relaxation time with temperature. 



The rate of electronic deexcitation from an upper state to a 
lower state (by electron impact in a molecule) is given by 

‘'nm " [ACT/10, 000)’^gj^]/Cm5 n5) [(1/P) . (l/P)] (95) 

where ff and IT are the principal quantum numbers of the lower and 
upper states, respectively, and A is the excitation rate constant 
(different for different molecules) . Due to the absence of e close 
to the wall, where the NLTE effect is more pronounced, deexcitation 
by electronic collisions is not considered. Here, only the colli- 
sional deactivation of by heavy particles is considered. 

The collisional relaxation time for a hydrogen ion is given 
(ref. 61) by the relation: 

n^(H ) = [m. (3kT) ]/(17.94 n^e^^Z^ InA) (96) 


where m^ represents mass of the ions, n^ is the number density 
of ions, and A is a parameter which is expressed as a function of 
temperature by 


A 


(k^T^/ Trn^) 


1/2 


(97) 


For hydrogen ions, equation (96) may be simplified further as 
n^(H'^) = (22.8 T3/2)/(n. InA) 


(98) 


where T is the heavy particle temperature in K. The collisional 
relaxation times for the hydrogen ions are shown in figure 6 for 
different number densities and temperature values. 
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THERMODYNAMIC AND TRANSPORT PROPERTIES 


Thermodynamic properties for specific heat, enthalpy, and free 
energy, and transport properties for viscosity and thermal conduc- 
tivity are required for each species considered in the shock layer. 
The general expressions for total enthalpy, specific enthalpy, and 
specific heat at constant pressure are given respectively by 

= h + + v^)/2 

h = IC. h. 

1 1 


C 

P 


zc. c . 

1 pi 


In the present study, values for the thermodynamic and transport 
properties are obtained by using polynomial curve fits. The expressions 
for h. and C . are given (refs. 62, 63) by 

1 pi ^ 


h. = RT [A +(B/2)T + (C/5)T^ + (D/4)t3 + (E/5)t'' 

+(F/6)T^ + (G/7)tS] 

C . = R [A + BT + CT^ + DT^ + ET^^ + FT^ + GT^] 

pi 

where R is the universal gas constant and T is the equilibrium 
fluid temperature in the shock layer. For different species i 
under the present investigation, the polynomial coefficients A, 

B, . . . G are given in tables 2 and 3. Equations (99) to (101) 
are used to calculate the enthalpy variation in the shock layer. 


(99) 

( 100 ) 
( 101 ) 


( 102 ) 

(103) 
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Table 2. Coefficients for evaluation of the specific heat at 
constant pressure for hydrogen/ helium species. 





Coefficients 




Species* 

A 

B 

C 

D 

E 

F 

G 

H 

.2500E+01^ 

.2500E+01^ 

.2475E+013 

.7366E-04 

-.2537E-07 

.2386E-11 

-.4551E-16 

.2547E+05 

.2547E+05 

.2523E+05 

-.4601E+00 

-.4601E+00 

-.3749E+00 

H 

2 

.3057E+01^ 

.3100E+01^ 

.3363E+01^ 

.2676E-02 
.5111 £-03 
.4656E-03 

-.5809E-05 

.5264E-07 

-.5127E-07 

.5521E-08 

-.3491E-10 

.2802E-11 

-.1812E-11 

.3694E-14 

-.4905E-16 

-.9889E+03 

-.8773E+03 

-.1018E+04 

-.2299E+01 

-.1962E+01 

-.3716E+01 

h'^ 

.2500E+01^ 

.2500E+01^ 

.2600E+01^ 





.1840E+06 

.1840E+06 

.1840E+06 

-.1153E+01 

-.1153E+01 

-.1153E+01 

HE 

.2500E+0r 

•2500E+01^ 

.2500E+01^ 

0.0 

0.0 

0.0 

0.0 

-.7453E+03 

-.7453E+03 

-.7453E+03 

.9153E+00 

•9153E+00 

.9153EtOO 

HE^ 

.2500E+01^ 

.2500E+01^ 

.2500E+01' 





.2853E+06 

.2853E+06 

.2853E+06 

.1608E+01 

.1608E+01 

.1608E+01 

E' 

.2500E+01^ 

.2500E+01'' 

.2508E+01^ 

-.-332E-05 

.1364E-08 

-.1094E-12 

.2934E-17 

-.7453E+03 

-.7453E+03 

-.7450E+03 

-.1173E+02 

-.1173E+02 

-.1208E102 

★ 

Each species has 

three sets 

of values val id 

over three ( 

different temperature range; 

these are 

1 

> 300 K 2 > 

1000 K 

3 > 6000 K 






an 

on 



U1 

Table 3. Coefficients for evaluating the specific heat at 
constant pressure for the ablative products. 


Coefficients 


Species 

A 

B 

C 

D 

E 

F 

G 

C H 
<♦ 

.4968E+01 

.6531E+0! 

.6531E+01 

.1727E-01 

.6506E-02 

.6506E-02 

-.2994E-04 

-.2251E-05 

-.2251E-05 

.3246E-07 

.3329E-09 

.3329E-09 

-.1366E-10 

-.1721E-13 

-.1721E-13 

.7545E+05 

.7535E+05 

.7535E+05 

-.8769E+00 

-.7446E+01 

-.7446E+01 

0 

.2946E+01 

.2542E+01 

.2546E+01 

-.1638E-02 

-.2755E-04 

-.5952E-04 

.2421E-05 

-.3102E-08 

.2471E-07 

-.1602E-05 

.4551E-11 

-.2798E-11 

-.1602E-12 

-.4368E-15 

.9380E-16 

.3890E-05 

.2923E+05 

.2915E+05 

.2963E+01 

.4920E+01 

.5049E+01 

0 

2 

.3526E+01 

.3621E+01 

.3721E+01 

-.1878E-02 

.7361E-03 

•4254E-03 

.7055E-05 

-.1965E-06 

-.2835E-07 

-.6764E-03 

.3620E-10 

.6050E-12 

.2155E-11 

-.2894E-14 

-.5186E-17 

-.1047E+04 

-.1201E+04 

-.1044E+04 

.4305E+01 

.3615E+01 

.3254E+01 

o" 

.2498E+01 

.2506E+01 

.2944E+01 

.1141E-04 

-.1446E-04 

-.4108E-03 

-.2976E-07 

.1244E-07 

•9156E-07 

.3224E-10 

-.4685E-11 

-.5848E-11 

-.1237E-13 

.6554E-15 

.1190E-15 

.1879E+06 
. 1879E+06 
.1879Et06 

.4386E+01 

.4347E+01 

.1750E+01 

CO 

.3790E+01 

•2984E+01 

.3366E+01 

-.1619E-02 

.1489E-02 

.8027E-03 

.3692E-05 

-.5789E-06 

-.1968E-06 

-.2031E-08 

.1036E-09 

.1940E-10 

.2395E-12 

-.6935E-14 

-.5549E-15 

-.1435E+05 

-.1424E+05 

-.1434E+05 

.2955E+01 

.6347E+01 

.4263E+0] 

CO 

2 

.2400E+01 

.4460E+01 

.4413E+01 

.8735E-02 

.3098E-02 

.3192E-02 

-.6607E-05 

-.1239E-05 

-.1298E-05 

.2002E-08 

.2274E-09 

.2415E-09 

.6327E-15 

-.1552E-13 

-.1674E-13 

-.4837E+05 

-.4986E+05 

-.4894E+05 

.9695E+01 

-.9863E-I-00 

-.7288E-100 



Table 3. Cconcluded) . 





Coefficients 




Species 

A 

B 

C 

D 

E 

F 

G 

C 

.2532E+01 

.2581E+01 

.2141E+01 

-.1588E-03 

-.1469E-03 

.3219E-03 

.3068E-06 

.7438E-07 

-.5498E-07 

-.2677E-09 

-.7948E-11 

.3604E-11 

.8748E-13 

.5890E-16 

-.5564E-16 

.8524E+05 
.8521 E+05 
.8542E+05 

.4606E+01 

.4312E+01 

.6874E+01 

C 

2 

.7451E+01 

.4043E+01 

.4026E+01 

-.1014E-01 

.2057E-03 

.4857E-03 

.8587E-05 

.1090E-06 

-.7026E-07 

.8732E-09 

-.3642E-10 

.4666E-11 

-.2442E-11 

.3412E-14 

-.1142E-15 

.9891 E+05 
.9970E+05 
.9787E+05 

-.1584E+02 

.1277E+01 

.1090E+01 

C 

3 

.5564E+01 

.4394E+01 

.2213E+02 

-.4687E>02 

.2964E-02 

-.1759E-01 

.1533E'*04 

-.1541E-05 

.5565E-05 

-.1509E-07 

.3796E-09 

-.6758E-09 

.5100E-11 

-.3010E-13 

.2925E-13 

.9926E+05 

.9926E+05 

.9423E+05 

-.2000E+01 

.2970E+01 

-.1021E+03 


.2593E+01 

.2511E+01 

.2528E+01 

--.4068E-03 

-.1735E-04 

.4869E-05 

.6892E-06 

.9504E-08 

-.7026E-08 

-.5266E-09 

-.2218E-11 

,1134E-11 

.1508E-12 

.1862E-15 

-.3476E-16 

.2166E+06 

.2166E+06 

•2168E+06 

.3895E+01 

.4286E+01 

.4139E+01 

C H 
2 

.2649E+01 

.4420E+01 

.5307E+01 

.8491E-02 

.22HE-02 

.8966E-03 

-.9816E-05 

-.5929E-06 

-.1378E-06 

.6537E-08 

.9419E-10 

.9251E-11 

-.1735E-11 

-.6852E-14 

-.2278E-15 

.5627E+05 

.5583E+05 

.5809E+05 

.7689E+01 

-.1158E+01 

-.5288E+01 

C H 
2 2 

.1410E+01 

.4575E+01 

..6789E+01 

.1905E-01 

.5123E-02 

.1503E-02 

-.2450E-04 

-.1745E-05 

-.2295E-06 

.1539E-07 

.2867E-09 

.1534E-10 

-.4134E-11 

-.1795E-13 

-.3763E-15 

.2618E+05 

.2560E+05 

.2590E+05 

.1139E+02 

-.3573E+01 

-.1539E+02 

C H 

3 

.3344E+01 

.3877E+01 

.3877E+01 

.1068E-01 

.6724E-02 

.6724E-02 

-.1331E--04 

-.2605E-05 

-.2605E'05 

.1338E-07 

.4416E-09 

.4416E-09 

-.5698E-11 

-.2708E-13 

-.2708E-13 

.6258E+05 

.6256E+05 

.6256E+05 

.6000E+01 

•3826E+01 

.3826E+01 

1 

> 300 K 2 > 

1000 K 3 

> 6000 K 






Cn 
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For the shock- layer gas; the mixture viscosity and thermal 
conductivity are obtained by using the semi-empirical formulas of 
Wilke (ref. 64) as; 


N 


i=l 1=1 ^ J 


IN IN 

K = E [XjK /(i; xj. )] 

1=1 ^ ^ j=l ^ 


where 


... = [1 + CupupV2 V*]^/ j/8 [1 + 


(104) 


(105) 


(106) 


and is the molecular weight of species i. For hydrogen/helium 

species, specific relations for viscosity and thermal conductivity 
are given in reference 65. The viscosity of and He, as a function 
of temperature, can be obtained from reference 35 as 


= (0.66 X 10"^) (T) 3/2/(t + 70.5), N sec m ^ 

Ho 


(107) 


= (1.55 X 10"^) (T) V2/(t + 97.8), N sec m ^ 
He 


(108) 


The specific relations for viscosity of other species are given in 
reference 35. The general relation for the thermal conductivity, 
based on equation (105) , is given (ref. 66) as 

K. = A + BT 
1 


( 109 ) 


The coefficients A and B for different species used in this 
study are given in table 4. 
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Table 4. Coefficients for thermal conductivity in equation (109"). 


Species 

COEFFICIENTS 

A 

B 

H 

.2496E-04 

.5129E-07 

H 

.3211E-04 

.5344E-07 

h" 

.2600E-03 


HE 

.2038E-04 

.3249E-07 

HE-^ 

.2600E-03 


E' 

.2600E-03 


C 

.2506E-04 

.7479E-08 

C 

.8590E-05 

.6233E-08 

C 

.6300E-05 

.5804E-08 

C+ 

.2600E-03 


C H 

.H26E-04 

.7439E-08 

C H 

.1126E-04 

.7439E-08 

C H 

,6300E-05 

.5804E-08 

C h 

.6300E-05 

.5804E-08 

0 

.1250E-04 

.7092E-08 

0 

.1019E-04 

.490iE-08 

0+ 

.2600E-03 


CO 

.8590E-05 

.6233E-08 

CO 

.8590E-05 

.6233E-08 



The heat transfer to the wall due to conduction and diffusion is 


referred to here as the convective heat flux and is given in terms 
of the thermal conductivity and viscosity as 

Ns 

q = -£2[K(9T/3n) + (jjLe/Pr) 12 OC./3n)h.] 

^ i=l ^ ^ 

The values for the Prandtl number Pr and the Lewis number Le are 
taken as 0.74 and 1.1, respectively. It should be noted that 
equation (110) is similar to equation (23). 


(110) 
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METHOD OF SOLUTION 


The numerical procedures for solving the inviscid and viscous 
flow equations are discussed in detail in references 5 and 29. 
Tiwari and Szema applied the method outlined in reference 5 in 
their study of the effects of precursor heating on chemical and 
radiative nonequilibrium viscous flow around a Jovian entry body 
(ref. 8) . A modified form of this procedure is used in this study 
to obtain solutions of the shock- layer equations under both the 
equilibrium and nonequilibrium radiative heat transfer conditions. 
In this method, a transformation is applied to the viscous shock- 
layer equations in order to simplify the numerical computations. 

In this transformation most of the variables are normalized with 
their local shock values. The transformed variables are 


n = n/n^ 

p = P/P3 

P = P/Ps 


? = s 

p = P/P3 

r = k/k 

s 


u = u/u 

s 

T = T/T 

s 

C = C /C 
P P Ps 


V = v/v 

s 

H = H/H 

s 


(111) 


The transformation relating the differentials are 

3/8s ( ) = 8/9? (l/n^)(dn^/d?)n 9/9n ( ) (112a) 


and 

9/9n ( ) = 1/n 9/9 ( ), 92 / 9^2 = i/n2 9^/9n2 ( ) (112b) 

ST) ^ 

The transformed equation can be expressed in a general form as 
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(3^W/3ri^) + (ai9W/8n)-+ Ca- 2 ''^) + (^ 3 ) + (a43W/9^) - 0 


(113) 


The quantity W represents u in the x-momentum equation, T 

in the temperature energy equation, H in the enthalpy energy equation, 

and C in the species continuity equations. The coefficients a^ 
i 

to a 4 are given in reference 5. With radiation included in the 
study, the coefficient is different from the expression given 

in reference 5. The modified value for the a 3 in the enthalpy 
energy equation is 


Pr 


as = 


p n^ 
r s 


P pH 
s s 




cos 9 


1 + 


n rp< 
s 


r + 


n cos 
s 


+ 


n PrP V V 
s s s 


e^p pH 
s s 



n PrPr 
s s 


£“^4 H p 
s s 



/ K 

^ I 1 + n nx 
\ s 


cos 0 


r + n cos 
s 



(114) 


where 




ps 


n Pr 
s s 




N 

p\ ^ 

i=l 


9C. 


u^pu 


1 \ on 


Pr 


fPr Pr 
s 


- 1 ) 


9p 

9n 


] 


2 - 2 

p p icpu 

s s 

1 + n riK 
s 


(115) 


other transformed equations are the same as given in reference 5. 

The surface boundary conditions in terms of transferred variables 

are 


p = 0, V = 0, T = T^ 


(116) 
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The transformed shock conditions are found to be 


u = v = T= H = p= p = l 
at n = 1 . 

The second-order partial differential equations as expressed by 
equation (113), along with the surface boundary and shock conditions, 
are solved by employing an implicit finite-difference scheme. The 
procedure is discussed briefly in the Appendix where flow diagrams 
for specific computations are also provided. 


(117) 
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RESULTS AND DISCUSSION 


Introduction 

The results obtained in the present study extend over a wide 
range of free-stream and flow conditions and probe shapes. The 
basic shapes considered for the entry probes are sphere cones, 
hyperboloids and ellipsoids. Most results have been obtained 
for the entry conditions (i.e., the free-stream conditions at dif- 
ferent entry altitudes) given in table 1. However, some specific 
results have also been obtained for other entry conditions. The 
shock- layer gas has been assumed to be in chemical equilibrium for 
the entire study. Different spectral models for radiative transfer 
in the gas have been considered, and various results obtained by 
these models are compared. The NLTE analysis was first carried out 
by considering only the hydrogen/helium species in the shock- layer 
gas. Later, the contributions of ablative products were also included 
in the NLTE analysis. Thus, in accordance with the four areas of 
this study, the results are presented in the following four subsections: 

(1) significance of radiation models on the flow-field solutions, 

(2) influence of NLTE radiation without ablation injection, [3) 
importance of NLTE radiation with ablation injection, and (4) effect 
of probe shape change on the flow phenomena. The physical model 
and flow conditions for which the results were obtained are given 

in each subsection. 

Significance of Radiation Models on the Flow-Field Solutions 

By employing the three different absorption models discussed 
under "Radiative Transport Models" [subsection titled "Spectral 
Model for Gaseous Absorption") results were obtained for the flow- 
field variables and the wall radiative heat flux distribution for 
different entry conditions and body configurations. Inviscid as 
well as viscous results were obtained for a 55-degree half angle 
sphere cone, while only viscous results were obtained for a 50- 
degree hyperboloid. Comparison of inviscid and viscous results 
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is presented for a 55-degree sphere cone entering the Jovian atmosphere 
at an altitude of Z = 116 km. Next, a series of viscous results is 
presented for the 55-degree sphere cone at different entry conditions. 
Finally, viscous results for a 50-degree hyperboloid are presented. 

Comparison of inviscid and viscous results.- Inviscid and viscous 
results obtained by employing the detailed and 30-step radiation 
models are compared in this subsection. The temperature distribution 
along the stagnation streamline is illustrated in figure 7(a). The 
agreement between inviscid and viscous results is seen to be fairly 
good except near the body, where viscous boundary- layer effects are 
predominant [see fig. 7(b)]. The difference between the detailed 
and 30-step model results is lower for the inviscid case than the 
viscous case. This is due to relatively higher temperature across 
the shock-layer for the inviscid analysis. As pointed out earlier, 
the step model is more accurate at higher temperatures. 

The shock standoff distance as a function of body location is 
illustrated in figure 8. The first three curves illustrate the 
inviscid results for the three different radiation models. The 
fourth curve, obtained by employing Nicolet’s detailed radiation 
model, is for the viscous case, and is drawn here for comparison. 

The shock standoff distance is slightly larger for the present 30- 
step model as compared with the results of the detailed and Sutton's 
58-step models. Although the difference between inviscid and 
viscous results is seen to be quite small, use of the viscous 
analysis is recommended for more realistic and accurate calculations. 

Results of radiative heating along the body are illustrated in 

figure 9. While inviscid results are seen to be slightly higher 

at the stagnation point, viscous results are relatively higher at 

* * 

other body locations (up to s = 0.6). This is a direct conse- 
quence of viscous boundary- layer effects. A discussion of viscous 
results for different radiation models follows. 



TEMPERATURE 



(a) /I = 

Figure 7. Temperature distribution along the stagnation streamline for inviscid and viscous 
analysis (55° sphere cone, Z = 116 km). (Continued) 










Viscous results for a 55- degree sphere cone. - Viscous results 
for a 55-degree cone (with a nose radius of 22.2 cm) for different 
entry conditions are now presented. Results of various radiation 
models are compared in order to establish the validity of the present 
30-step radiation model. 

In the shock layer, the temperature distribution along the 
stagnation streamline is illustrated in figure 10 for two different 
free-stream (density) conditions. It is found that the present 30- 
step model underpredicts the shock- layer temperature by a maximum 
of 11 percent in comparison to Nicolet's detailed model and by about 
4.5 percent when compared with Sutton's 58-step model. For free- 
stream conditions resulting in higher shock temperature, the agreement 
between the results is even better. This is because the higher 
temperature absorption spectrum can be approximated accurately 
by the present step-model. 

From the results presented in figures 7 and 10, it is noted 
that there exists a steep temperature gradient in the regions close 
to the body. At locations about five times the nose radius (normal 
to the body), only a slight variation in the shock-layer temperature 
is noticed. This fact was utilized in dividing the shock-layer into 
different temperature zones for evaluating the absorption coefficient. 
In a preliminary study, two methods were used to account for the 
temperature dependence of the absorption coefficient. In the 
first method, the absorption coefficient was calculated at the 
shock temperature T*. This value was used in analyzing the flow 
field in the entire shock layer. Results obtained by this method 
are designated here as "present-approximate" results. In the second 
method, the shock- layer is divided into three different temperature 
zones, two of which are closer to the body (because of the steep 
gradient near the body). For each temperature zone, a different 
30-step model for absorption is obtained. These are read as input 
in the computer program while evaluating the flow variables in the 
particular temperature zone. Results obtained by this method are 
denoted here as "present" results. 
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The variation in temperature just behind the shock (at location 
n = 0.05) with distance along the body surface is illustrated in figure 
11 for entry conditions at Z = 116 km. The results of the present 
model are found to be about six percent lower than the results of 
Nicolet's model. This difference is seen to be fairly uniform along 
the body. 

Figure 12 shows the shock standoff variation with distance along 
the body surface for entry conditions at Z = 116 km. Results of 
Sutton's model are found to be in general agreement with the results 
of Nicolet's model. The present model is seen to overestimate the 
results by a maximum of 8.6 percent when compared with the results 
of Nicolet's model. This is mainly because the present model 
underpredicts the shock-layer density. 

The radiative heating rate along the body surface is illustrated 
in figures 13 and 14 for different entry conditions. As would be 
expected, in all cases, the maximum heating occurs at the stagnation 
point. For Z = 116 km, results presented in figure 13 show that the 
present model underpredicts the heating rate by a maximum of 13.6 
percent when compared with Nicolet's model. For the case of higher 
free-stream density (and hence a higher shock temperature) , differences 
in the results of the present and other models are seen to be smaller. 
Figure 14 shows the results of radiative heating for 131-km entry 
conditions. For this higher altitude, the heat transferred to the 
body is lower because of lower free-stream density and pressure. 

For this case, differences in the results of the present and Nicolet's 
model are seen to be slightly higher. 

Viscous results for a 50-degree hyperboloid. - Viscous results 
for a 50-degree hyperboloid (with a nose radius of 22.2 cm) are 
presented in this subsection for different entry conditions. The 
temperature distribution in the shock layer (along the stagnation 
streamline) is illustrated in figure 15 for entry conditions at 
Z = 116 km. The results of the three radiation models are seen 
to follow the same general trend as for the 55-degree sphere cone. 

A maximum difference of about 4.5 percent is seen between the 
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Figure 11. Temperature variation just behind the shock (n = 0.05) with distance along the body 
surface (55° sphere cone, Z = 116 km). 
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Figure 12. Shock standoff variation with distance along the body surface (55° sphere cone 
Z = 116 km). 
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Figure 13. Radiative heating along the body for two different free-stream 
densities (55® sphere cone) . 
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Figure 14. Radiative heating along the body for entry conditions at Z = 131 km (55*^ sphere cone) 
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Figure 15. Temperature distribution along the stagnation streamline 
(50° hyperboloid, Z = 116 km). 


present model and Nicolet's jnodel. This difference is near the 
body (at n = 0.0095). Agreements between the results are better 
towards the shock. 

The variation in temperature just behind the shock (at location 
n = 0.07) with distance along the body surface is illustrated in 
figure 16 for entry conditions at Z = 116 km. The results indicate 
very good agreement between the three radiation models. The results 
of the present model are within 1.4 percent of the results of 
Nicolet's model. As would be expected, maximum difference in results 
occurs at the stagnation streamline. 

The shock standoff variation with distance along the body surface 
for entry conditions at Z = 116 km is shown in figure 17. As was the 
case with the 55-degree sphere cone, the present model is seen to 
overestimate the results in comparison to the other models. 

The radiative heating rate along the body surface is illustrated 
in figures 18 to 20 for different entry conditions. For this body 
geometry also, the maximum heating occurs at the stagnation point. 

For Z = 116 km, results presented in figure 18 indicate that the 
present model underpredicts the heating rate by a maximum of about 
13 percent when compared with Nicolet's model. For higher free-stream 
density, the results presented in figure 19 show smaller differences 
in the results of various radiation models. For entry conditions 
at Z = 131 km, results presented in figure 20 indicate that heat 
transferred to the body is significantly lower. This is because of 
lower free-stream density and pressure. As was the case with the 
55-degree sphere cone at this altitude, the difference between the 
present and Nicolet's results is relatively higher. 

It is found that use of the present model reduces the computa- 
tional time significantly. The use of this model is recommended for 
simple parametric study. However, the use of the present model has 
its limitations. For problems with varying shock-layer compositions 
and large temperature variations, a suitable empirical correlation 
has to be developed from the detailed model to make the present model 
more versatile and accurate. 
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Figure 17. Shock standoff variation with distance along the body surface (50° hyperboloid, Z = 
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Figure 19. Radiative heating along the body for £ree-st: 
(50° hyperboloid) . 
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Influence of NLTE Radiation Without Ablation Injection 


Specific results obtained by using the NLTE formulation for 
the radiative transport, as described by equations (64) to (66), 
are presented in this subsection. The governing equations (6) to 
(13) were used to solve the shock- layer flow which was considered 
to be viscous and in chemical equilibrium. The effect of ablation 
mass injection into the flow was neglected. In this subsection, 
equilibrium radiative transport solutions are presented along with 
the NLTE results for comparative purpose. For this study, the entry 
bodies considered are a 50-degree hyperboloid and a 55-degree sphere 
cone which enter the Jovian atmosphere at a zero-degree angle of 
attack. In both cases, the body nose radius is taken to be 

22.2 cm. The body surface is assumed to be gray having a surface 
emittance of 0.8, and the wall temperature T^ is taken to be 
constant at 4,200 K. The variation of the nonequilibrium absorption 
coefficient (as compared to the equilibrium values) is shown in 
figure 21 at a temperature of 15,950 K. Within the confines of 
assumptions made in this study, the equilibrium and nonequilibrium 
absorption coefficients are found to be about the same beyond 10 eV. 
The shaded portion in the figure represents the decrease in the 
absorption coefficient values as a result of an increase in the 
population of higher energy levels. 

NLTE results for a 50-degree hyperboloid . - The variation of 
collisional relaxation time across the shock layer is shown in figures 
22(a) to 22(c) for three different entry conditions. Figure 22(a) 
shows the results for combined H-H and 0.5 H - 0.5 collisions, 
figure 22(b) for H-H and H^-H^ collisions, and figure 22(c) for 
0.95 H - 0.05 H^ and H^-H^ collisions. The results indicate that 
in all three cases the relaxation time does not vary significantly 
across the shock layer except very close to the body. Thus, higher 
NLTE effects will be expected in regions closer to the wall where 
the assumption of chemical equilibrium usually is justified (ref. 8) . 
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Figure 22. (Continued]. 
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Figure 22, (Concluded). 
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Since relaxation times are comparatively longer for higher entry 
altitudes, one would expect the NLTE effects to be greater at higher 
altitudes. The results presented in figure 22(a) for combined H-H 
and 0.5 H - 0.5 collisions, of course, are not representative of 
the results for the actual shock-layer collisional processes; they 
are presented here only for comparative purposes. The results of 
the other two collisional processes are found to be very close 
[figs. 22(b) and (c)], and they do represent very nearly the results 
of actual shock- layer collisional processes. 

For the entry condition at Z = 116 km, the temperature distri- 
bution along the stagnation streamline is illustrated in figures 23(a) 
and (b) . The results in figure 23(a) are for combined H-H and H^-H^ 
collisions and in figure 22(b) for 0.95 H - 0.05 H 2 and H^-H^ colli- 
sions; the two sets of results are found to be almost identical. 

The results clearly indicate that the NLTE temperature distribution 
is consistently lower than the equilibrium temperature. This implies 
that the shock-layer gas absorbs less energy under NLTE conditions 
than under equilibrium conditions. This is because under NLTE 
conditions (where the population ratios of the energy levels deviate 
from the equilibrium Boltzmann distribution) the number of particles 
(capable of absorbing the incoming radiation) in the ground state 
is comparatively less than under equilibrium conditions. The dif- 
ferences between LTE and NLTE results are seen to be lower toward 
the shock than toward the body. This is mainly due to direct dependence 
of the collisional deactivation process on the temperature. The 
maximum NLTE effect, therefore, will occur near the body surface 
where the collisional deactivation process is slower because of 
the lower temperature (see also the results presented in fig. 22). 

For the entry conditions at Z = 116 km, the temperature distribu- 
tions along the body are shown in figure 24 for two locations in the 

shock layer. The NLTE results were obtained by considering the combined 
+ + , ^ 

H-H and H -H collisional process. In the region close to the shock 
(n = 0.068), the differences between the LTE and NLTE temperatures are 
negligible (less than 0.1 percent). On the other hand, at a location 
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Figure 24. Temperature variation along the body for tv;o locations in 
the shock layer, H-H and collisions. 


closer to the body (n = 0.00012), the differences are significantly 
larger: a difference of about 6.8 percent is noted at the stagnation 

point. Since temperature decreases in the direction of the s coor- 
dinate, the NLTE influence is found to be greater at locations away 
from the stagnation streamline. 

For the entry conditions at Z = 116 km, the shock-standoff variation 
with distance along the body surface is shown in figure 25. Since the 
shock-standoff distance is influenced by the entry conditions and the 
shape of the entry body, the conditions of NLTE in the shock layer do 
not have any effect on its variation. 

For entry conditions at Z = 109 km, the LTE and NLTE radiative 
heating along the body are illustrated in figure 26. The NLTE results 
were obtained by considering the combined H-H and collisional 

process. The results simply indicate that the NLTE heating is con- 
sistently lower than the LTE heating all along the body. Since the 
number density of participating particles is relatively higher at lower 
altitudes, larger NLTE effects would be expected at altitudes higher 
than Z = 109 km. 

For the peak heating entry conditions (i.e., for Z = 116 km), the 
LTE and NLTE results of radiative heating along the body are illustrated 
in figure 27. In order to assess the influence of various deactivation 
processes, the NLTE results have been obtained by considering five 
different collisional relaxation times. This is essential because the 
exact nature of the collisional deactivation process, which actually 
occurs in the shock-heated gas, is not known. It is evident from the 
figure that the NLTE results obtained by considering only the 
collisions are very close to the LTE results. Consequently it may be 
concluded that in a fully ionized plasma the assumption of LTE is 
justified. The NLTE results obtained by considering only the 0.5 H 
- 0.5 H 2 collisions are seen to be significantly lower than the LTE 
results. This, however, does not represent a physically realistic 
situation for the shock-layer gas (because of a very low number density 
of hydrogen molecules); the results are presented here only for 
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Figure 26. Radiative^heating along the body for entry conditions at Z = 109 km, 
H-H and H -H collisions. 
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Figure 27. Radiative heating along the body for entry conditions at Z = 116 km. 
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comparative purposes. Probably the more realistic collisional process 
for the shock-layer gas may be represented by the combined H-H and 

or H-H 2 collisions. The NLTE results obtained by considering the 
combined relaxation times of these collisions are seen to be much lower 
than the LTE results all along the body. The maximum NLTE effects are 
found to be for the case of combined H-H and 0.95 H - 0.05 H 2 collisions. 
The difference between LTE and NLTE results for this case is 11 percent 
at the stagnation point. The case of the combined H-H and H^-H^ col- 
lisional process (which is very close to the case of the combined 
0.95 H - 0.05 H 2 and H^-H^ process), however, appears to be more pliys- 
ically convincing. For this case, a comparison of results presented in 
figures 26 and 27 reveals that NLTE effects are significantly higher 
at Z = 116 km than at Z = 109 km. Figure 27 shows a nine percent 
reduction in radiative heating for this case at the stagnation point. 
Perhaps an even more convincing process to consider will be the combined 
H-H, H -H , and 0.95 H - 0.05 H 2 collisional process. However, at 
present no theory is available to calculate the relaxation times of such 
collisions. The effects of such collisions, of course, will be lower 
than those for the combined H-H and 0.95 H - 0.05 H 2 collisions; the 
results are expected to be closer to the results of the combined H-H 
and H^-H^ collisions. 

For entry conditions at Z = 131 km, the results of radiative and 
convective heating along the body are illustrated in figure 28. The 

+ 4 “ 

NLTE results were obtained by considering the combined H-H and H -H 
collisional process. For radiative heating, NLTE results are shown 
also for two-level energy transitions. The results clearly indicate 
that, although differences between LTE and NLTE results are small 
for the convective heating, they are considerably larger for the 
radiative heating. The contributions of higher level energy trans- 
itions on NLTE results are seen to be quite small (less than 1.7 
percent) . Since NLTE affects the convective heating only through a 
different temperature distribution, the effects are seen to be quite 
small away from the stagnation point. 
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Figure 28. Radiative and convective heating along the body for entry conditions at Z - 131 km, 
H-H and collisions. 
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Results of LTE and NLTE radiative heating at the stagnation point 

are shown in figure 29 for entry conditions at different altitudes. 

The NLTE results again are obtained by considering the combined H-H 
^ ■ ■ 

and H -H collisions. Once again, the NLTE results for two-level 
energy transitions are presented for comparative purposes. The results 
show that differences between LTE and NLTE heating rates are larger at 
higher altitudes. As mentioned earlier, this is because the densities 
of participating species are relatively lower at higher altitudes (lower 
pressures), and this, in turn, results in longer collisional deactivation 
times . 

NLTE results for a 55-degree sphere cone . - For the 55-degree 
sphere cone, the radiative heating results for the peak heating condi- 
tions are illustrated in figure 30. The NLTE results were obtained by 
considering the combined H-H and H -H collisional process. Because of 
the numerical instability, it was possible to obtain accurate NLTE 
results only up to the tangency point. Obviously, further work is 
needed to improve the numerical procedure for obtaining the NLTE results 
toward the downstream regions. In the stagnation region, the results 
for the sphere cone show essentially the same trend as for the hyper- 
boloid. At the stagnation point, the difference between the LTE and 
NLTE results for the sphere cone is about the same as for the hyper- 
boloid. This clearly indicates that the NLTE results are not 
influenced significantly by the changes in forebody shapes. However, 
further work is needed to make this a definite recommendation because 
the LTE heating rates are influenced by the shape change (refs. 26, 27). 

The results presented herein indicate that, although the relaxation 
times for collisions between neutral particles decrease with increasing 
temperature, the reverse is true for the charged particles. It is also 
noted that the physically realistic collisional process for the shock- 
layer gas (in the absence of ablative products) is the combined H-H and 
H^-H^ deactivation process. Specific results indicate that NLTE effects 
are greater closer to the body than near the shock. This is because the 
NLTE results are influenced strongly by the temperature distribution in 
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the shock layer and the effects are lower at higher temepratures . It 
is further noted that the contribution of higher level energy transitions 
on the NLTE results is relatively small for all entry conditions. It 
is found that the influence of NLTE, in general, reduces the convective 
and radiative heating to the entry body. Although this effect is small 
for the convective heating, the radiative heating is influenced signif- 
icantly. The NLTE effects are greater for higher altitude entry 
conditions. A qualitative comparison of the results for a 50-degree 
hyperboloid and a 55-degree sphere cone (under identical physical and 
entry conditions) shows that the NLTE results are not influenced 
significantly by the change in the forebody configurations of the entry 
probe. 
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Importance of NLTE Radiation with Ablation Injection 


To investigate the importance of NLTE radiation, the results 
were obtained only for a 35-degree hyperboloid (with coupled ablation 
mass loss from a carbon-phenolic heat shield). This is because, 
for this case, LTE results were already available in the literature 
(refs. 36 - 39). For comparative purposes, selected NLTE results 
were also obtained for the case with no ablative products (i.e., 
for m = 0) . The NLTE results, in this case, were obtained by con- 
sidering the combined collisional deactivation process of H-H and 

+ -f* 

H -H . Viscous shock-layer results obtained for the peak-heating 
conditions are presented in this subsection. 

As discussed earlier (under "Radiative Lifetimes and Collisional 
Processes for the Shock-Layer Gases"), for NLTE study, it is essen- 
tial to know the nature of the collisional deexcitation processes 
and relaxation times of different shock- layer species in presence 
of the ablative products. Figure 31 illustrates the important 
species concentrations near the wall influencing the collisional 
process. In general, the molecules are concentrated near the 
wall and the number density rapidly reduces away from the wall as 
they dissociate into and atomic carbon. 

In the presence of the ablative products, the relaxation times 
for the combined collisional process of C^-C, C2-H, and H-H 

were used in the present study. The majority of NLTE results were 
obtained by using the combined radiative lifetimes of the Swan, 
Freymark, and Mulliken band system as indicated by (III) 
under "Radiative Lifetimes and Collisional Processes for the Shock- 
Layer Gases." However, some results for radiative heating rates 
were also obtained by using the other radiative lifetimes of the 
band system as indicated by (I), (II), and (IV). 

The temperature variation across the shock- layer (for loca- 
tion s = 0) is shown in figure 32 for both LTE and NLTE (III) con- 
ditions. Results with no mass injection are also shown here for 
comparison. As would be expected, the shock- layer temperature. 




Figure 31. Species concentrations in the vicinity of the wall. 
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n = n /R^ 


Figure 32. Temperature distribution across the shock layer. 



in general, is lower in the vicinity of the body in the presence 
of the ablative products (refs. 36 - 39). It is seen that the NLTE 
temperature distribution is lower than the equilibrium values through- 
out the shock layer. A maximum difference of 5.48 percent is noticed 
between the 2 values at n = 0.13. The C 2 molecules in the ablation 
layer (a region in the vicinity of the wall), which absorb less under 
NLTE conditions, possess less energy than the equilibrium value. This, 
in turn, results in lower temperature values in the ablation layer, 
and the trend continues in the entire shock layer. 

Figure 33 illustrates the density and enthalpy variations across 
the shock layer for LTE and NLTE (III) conditions. The enthalpy 
variation has a similar trend as the shock- layer temperature shown 
in figure 32. It was found that NLTE essentially had no influence 
on the pressure distribution in the shock layer. The density, 
however, is seen to be significantly higher for the NLTE case. 

This is a direct consequence of relatively lower NLTE temperatures 
in the shock layer. A maximum increase in density of about 5.5 per- 
cent is noticed at n = 0.15. 

The equilibrium and nonequilibrium shock-standoff variation 
with distance along the body surface is shown in figure 34 for cases 
with and without ablation injection. As noted earlier, the shock- 
standoff distance is not influenced significantly by the NLTE 
conditions for the case with no ablation injection. For the case 
with ablation injection, however, the NLTE (III) results are com- 
paratively higher than the LTE results. A possible reason for this 
behavior is the combination of enthalpy and density variation 
in the shock layer along with the energy loss at the shock for 
nonequilibrium conditions. 

Variations in the nondimensional surface pressure and heating 

rate along the forebody of the probe are illustrated in figure 35. 

These quantities are nondimensionalized by their respective stagnation 

values of p = 6.309 atm, q (LTE) = 201.849 MW/m^, and q 

^w,o ^w,o ^ ^ ^w,o 

(NLTE) = 208.927 NW/m^. It is seen that NLTE virtually has no 
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Figure 34. Shock standoff variation with distance along the body surface 





Figure 35. 


Pressure and heating rate variations along the body surface 



influence on the pressure distribution. However, the total heating 
rate (convective plus radiative) is increased significantly under 
NLTE conditions. The main contribution to the total surface heating 
was found to be the radiative heating. The explanation of increased 
NLTE radiative heating, in this case, is given in the discussion of 
results presented in figure 36. 

The results of radiative heating rates for different conditions 
are given in table 5 and are shown in figure 36. The LTE and NLTE 
(III) results are compared in figure 36(a) for the cases with and 
without ablation injection. The results clearly indicate that the 
radiative heating to the body, in general, is reduced significantly 
in the presence of the ablative products. For the case with no 
ablation, the NLTE results are found to be significantly lower than 
the LTE results; a decrease of about nine percent is noted at the 
stagnation point. In the presence of the ablative products, however, 
the results presented in figure 36 and table 5 show that NLTE results 
are comparatively higher than the LTE results. The reason for this 
is as follows: Under NLTE conditions, the number of molecules 

in the ground state (that are capable of absorbing the incoming 
radiation from the shock- layer gases) is less as compared to the LTE 
values (i.e., the number based on the Boltzmann distribution). This 
increases the transparency of the ablation layer which, in turn, 
results in higher heating of the entry body. This reverse trend 
in the NLTE heating rate is an important finding of this study. 

The results for the NLTE heating rate obtained by considering 
different radiative lifetimes are illustrated in figure 36(b). 

The results for cases (II) and (HI) were found to be 

about the same for all body locations (see table 5) . The results 
for n^ (I) are seen to give the smallest increase in NLTE heating 
whereas the results for (IV) provide the maximum heating rate 

to the body. These NLTE results, however, do not differ from each 
other considerably. The maximum increase in the stagnation-point 
heating is found to be about 3.5 percent for (HI) and about 

5 percent for (IV). Thus, based on the information of radiative 
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Figure 36(a). Variation of radiative heating rate along the body surface 
(with and without ablation) . 
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lifetimes of the band system available at the present time, it 
may be concluded that the NLTE effects will be maximum for the 
combined radiative lifetime represented by (IV), and this will 

increase the radiative heating to the body by a maximum of about 
five percent. 

The ablation mass loss rate from the body surface is shown 
in figure 37. As would be expected, the mass loss rate is higher 
under the NLTE conditions (because of the increased heating rates 
to the body) . 

In conclusion, the results presented here indicate that the 
temperature and enthalpy distribution in the shock layer are lower 
under the NLTE conditions. The NLTE increases the density in the 
shock layer, but it has no influence on the pressure variation. 

The radiative heating to the entry body is increased significantly 
because of NLTE and this, in turn, results in increased mass loss 
from the body. 


Table 5. Wall radiative heat flux for the case with 
ablation under LTE and NLTE conditions. 




Wall Radiative 

Heat Flux 

q^, MlV/m^ 


s = S*/Rf5 

q*(LTE) 

* 

q^CD 

q^CII) 

q^dii) 

q^Civ) 

o 

o 

201.849 

206.774 

209.017 

208.927 

211.987 

CM 

o 

188.172 

193.127 

196. 020 

195.891 

201.497 

0.4 

164.709 

169.658 

173.278 

173.127 

178.165 

0.6108 

135.040 

139.244 

141.497 

141.543 

144.782 

0.7854 

110.900 

114.765 

115.900 

115.974 

120.841 

1.200 

75.012 

78.767 

80.256 

79.012 

80.572 

1.500 

58.271 

61.813 

61.570 

61.902 

62.090 
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Figure 37. Ablation mass loss rate along the body surface. 
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Effect of Probe Shape Change on the Flow Phenomena 


The shock shape and shock- standoff distance are influenced 
greatly by the shape of the entry body. To study the influence 
of shape change on flow phenomena around the entry body, it is 
essential to specify several speculative profiles that will emerge 
from mass loss and recession of the heat shield just after the 
peak heating phase of the entry mission. This is important because 
the main data collection phase of the mission starts after the 
heating phase. At this stage, it is essential to consider the actual 
shape of the probe rather than the initial shape in investigating 
the flow field and aerodynamic stability of the entry body. It is 
necessary, therefore, to have analytic expressions that are capable 
of generating such desired shapes. A general relation for the 
shape change can be given by a quadratic form as 

y = 2Rj^x - Bx2 (118) 

where B represents the bluntness factor which determines the body 
shape. For negative values of B, the resulting shapes are a 
family of hyperbolas. For B = 1, circular or spherical shapes are 
obtained. A family of parabolas is obtained for B = 0, and positive 
values of B give elliptical shapes. 

Equation (118) is used in generating different shapes for the 
entry probe. The initial body shapes considered for this study 
are a 45-degree sphere cone (i.e., a 45-degree half-angle, spherically 
capped, conical body), a 35-degree hyperboloid (i.e., a hyperbolic 
forebody shape with an asymptotic angle of 9 = 35®), and a 45-degree 
ellipsoid. The reason for selecting a 35-degree hyperboloid (instead 
of a 45-degree hyperboloid) is that the mass losses for this and the 
45-degree sphere cone and 45-degree ellipsoid are comparable. For 
all initial shapes, the nose radius considered is 31.12 cm and the 
base radius is taken to be twice the nose radius. The final shape 
after the heating phase will depend upon the extent of absorption 
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of incoming radiation by the ablative products. If no radiation 
blockage is assumed, then it is possible to have severe recession 
of the forebody as well as of the afterbody. A 25 or 50 percent 
radiation blockage will result in relatively less severe blunt jj'ig 
of the entry probe. It is possible to generate these speculative 
profiles with the help of equation (118). 

The numerical procedure employed by Sutton et al . (ref. 26) 
for inviscid radiating flow is basically used in this study also. 

For the initial profile, the x and y coordinates, the distance 
along the body, and the radius of the curvature are specified at 
14 selected stations along the body. For the changed profiles, 
the distance along the body, the radius of curvature, and the body 
angle are calculated at each station by using a subroutine with 
X and y coordinates as inputs. A three-point central differencing 
scheme is used for calculation of the s location, and a two-point 
backward differencing scheme is used for calculation of the curvature. 

The three initial body shapes and the corresponding body profiles 
that emerge when the influence of shape change is considered are 
shown in figures 38(a) to (c) . Figure 38(a) shows the forebody 
configurations of the 45-degree sphere cone in which profile 1 
represents the initial shape and profile 2 is the corresponding 
blunted profile. Profile 2 represents a case where the entry body 
has experienced a severe blunting near the stagnation region as 
compared to the downstream region. Specifically, this represents 
a shape where the initial body has undergone a severe mass loss near 
the nose after absorbing about 50 percent of the incoming radiation. 
Profiles 3 and 4 in figure 38(a) represent two arbitrary blunted 
shapes for which the mass loss is assumed to be uniform all along 
the body. For all the forebody profiles shown in figure 38(a), 
the afterbody shape is a 45-degree half-angle cone. 

Different configurations for the 35-degree hyperboloid are shown 
in figure 38(b). In this case also, the severely blunted profile 
is represented by the second curve. This corresponds to the case 
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(c) for 45-degree ellipsoid. 


Figure 38. (Concluded). 


119 



of severe stagnation region mass loss with assumed radiation blockage 
of about 50 percent. Profile 3 in figure 38(b) represents the shape 
for 40 percent radiation blockage with uniform mass loss along the 
entire body. Different configurations for the 45-degree ellipsoid 
are shown in figure 38(c), and they correspond to exactly the same 
conditions as for the 35-degree hyperboloid. 

The condition of LTE for the radiative transport was assumed 
while obtaining the results using different probe configurations. 

Since NLTE is a condition of the absorbing/emitting gases, it is 
assumed that the shape change will not influence the NLTE phenomena. 
The results have been obtained for different entry conditions, and 
they are presented here first for the 45-degree sphere cone and 
then for the 35-degree hyperboloid. Finally, the peak heating 
results for the 45-degree sphere cone and 35-degree hyperboloid 
are compared with peak heating results of the 45-degree ellipsoid. 

Results for 45-degree sphere cone . - Different results obtained 
for the 45-degree sphere cone are illustrated in figures 39 to 46. 

The shock-standoff variation with distance along the body surface 
is shown in figures 39 and 40. For peak heating conditions (i.e. 
for Z = 116 km), results obtained for the four profiles indicated 
in figure 38(a) are illustrated in figure 39. Results obtained 
for the initial and blunted profiles are compared in figure 40 
for the tliree entry conditions considered. The results indicate 
that blunting of the nose region increases the shock- layer thickness 
not only near the stagnation region but all along the body, although 
there is no significant change in the body shape near the flank 
region. It is evident from figure 39 that profiles 3 and 4 do not 
influence the shock- standoff distance appreciably. This is because 
the uniform mass loss tends to preserve the original configuration 
of the entry body. As would be expected, the shock-standoff dis- 
tance, in general, increases with decreasing altitude, and near the 
stagnation region the influence of nose blunting is greater at 
lower altitudes. 
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Figure 39. Shock-standoff variation for 
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Figure 40, 


Shock-standoff variation for 45-degree sphere cone. 


Variation of the shock temperature (i.e., the temperature just 
behind the shock) with distance along the body surface is shown 
in figure 41 for the three entry conditions considered. For these 
conditions, the shock temperature, in general, is found to be higher 
for higher entry altitudes. The effect of blunting is seen to increase 
the shock temperature along most of the downstream region of the probe. 
For Z = 116 km, a maximum difference of 5.6 percent is found at the 
tangency point (where the forebody and the afterbody coincide asymptot- 
ically). The effect of blunting is seen to have relatively higher 
influence on the shock temperature for the other two entry conditions. 
Results of moderate shape change (with uniform mass loss) , as shown 
by curves 3 and 4, indicate negligible influence on the shock tempera- 
ture. Variation in density just behind the shock along the body sur- 
face is shown in figure 42 for 2 entry conditions, Z = 116 km and 
138 km. In general, lower shock densities are associated with higher 
entry altitudes. This is because the free-stream densities are lower 
at higher altitudes. The shock-density variation is relatively higher 
for the blunted profile, and a significant increase is noted from 
the stagnation point to the tangency point. However, virtually no 
difference is seen after the location s = 1. For Z = 116 km, a max- 
imum difference of 5.6 percent is noted between the initial and 
blunted profiles at location s = 0.8. This difference is even smaller 
(4.5 percent) for results at Z = 138 km. Results of profiles 3 and 
4 are in general agreement with the results of the initial profiles; 
a maximum difference of 1.25 percent is noted at s = 0.8 for Z = 116 km. 

Variations in density and v-velocity across the shock layer are 
shown in figures 43 and 44, respectively. Figure 43 shows the density 
variation for body locations (s = 0 and 1.4) and entry conditions 
(Z = 116 km and 138 km) . It is seen that along the stagnation line 
the density is not influenced by the shape change. However, signifi- 
cant differences in results of the initial and blunted profiles are 
noted for the downstream location of s = 1.4. The density values 
are lower for the blunted profile because the temperatures are relatively 
higher. Uniform profile changes do not alter the density values 
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appreciably. Figures 44(a) and (b) show that shape change does not 
influence the radial component of velocity significantly along the 
stagnation streamline; however, blunting is seen to influence the 
velocity toward downstream locations. 

In the shock layer, variation in pressure is relatively small 
as compared to variations in temperature and density. The variation 
of pressure along the wall is shown in figure 45 for the four profiles 
and three entry conditions. As would be expected, pressure distribu- 
tion is relatively higher for lower altitudes, and the maximum pressure 
occurs at the stagnation point. Results of profiles with uniform 
mass loss are not significantly different from the results of the 
initial profile. Blunting of the entry body is seen to increase the 
wall pressure significantly, the increase being maximum closer to the 
tangency point. Blunting does not seem to affect the stagnation 
region and downstream pressure distribution appreciably. However, 
it is possible for the blunted probe to experience relatively higher 
total drag. 

The radiative heating results are illustrated in figures 46(a) 
and (b) . In each case, the heating rate is seen to be significantly 
higher for the blunted profile all along the body. This, however, 
is expected because the shock temperature and shock-standoff distance 
are relatively higher for the blunted profile (see figs. 40 and 41). 

The maximum stagnation point heating occurs, of course, at Z = 116 km. 
For this case, the blunted profile heating rate is about nine percent 
higher than the initial profile heating. Figure 46(a) shows that 
heating rates for profiles 3 and 4 are lower than those of the initial 
profile. A difference of 4.9 percent between initial profile and 
profile 3, and of 5.3 percent between initial profile and profile 4, 
is noticed at the stagnation point. Results presented in figure 
46(b) indicate that the stagnation-region heating is comparatively 
higher for Z = 109 km than for Z = 138 km. The difference between 
stagnation-point heating rates for the blunted and initial profiles 
is 10 percent for Z = 109 km and 6 percent for Z = 138 km. The 
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Figure 45. Surface pressure variation for 45-degree sphere cone. 
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results clearly indicate that the shape change can have significant 
influence on heating of the afterbody of the entry probe. 

Results for 55-degree hyperboloid . - Some important results 
obtained for the 35-degree hyperboloid profiles [see fig. 38 [b)] 
are presented in figures 47 to 53. As expected, the hyperboloid 
results show a smoother trend than the sphere cone results. The 
shock-standoff variation with distance along the body surface is 
illustrated in figure 47 for the three entry conditions considered. 

The results indicate that shape change increases the shock-standoff 
distance all along the body. The increase is greater for the blunted 
profile, and a maximum increase of 8.5 percent is noted for Z = 109 km 
at s = 0. The results of profile 3 show an increase of only about 
3 percent for Z = 109 km at s = 0. Variation of the shock temperature 
is illustrated in figure 48. The results show that the shape change 
has only a slight influence on the temperature between locations 
s = 0.2 and 1.2, and its effects are negligible further downstream. 
Variation in the shock density for the initial and blunted profiles 
is shown in figure 49 for entry conditions at Z = 116 km and 138 km. 

The results indicate the shape change has only a slight influence 
on the density variation. 

The changes in velocity, density, temperature, and pressure 
across the shock layer of a hyperboloid essentially follow the same 
general pattern as for the sphere cone, but the effects of shape 
change are not as pronounced. Variations in density and v-velocity 
across the shock layer are shown in figures 50 and 51, respectively. 
Figure 50 shows the density variation for two body locations (s = 0 
and 1.5] and entry conditions (Z = 116 and 138 km). It is seen 
that along the stagnation line the density is not influenced by the 
shape change. However, small differences in results of the initial 
and blunted profiles are noted for the downstream location of s = 1.5. 
Uniform profile changes alter the density only slightly in the regions 
closer to the body. Figure 51 shows that along the stagnation streamline 
the shape change does not have any influence on the radial component 
of velocity: only slight changes are noted for Z = 138 km and s = 1.5. 
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Figure 48. Variation of temperature just behind the shock for 35-degree hyperboloid. 
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Figure 50. Variation of density in the shock layer for body locations s = 0 and 1.5, 
35-degree hyperboloid. 
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Figure 53. Radiative heating rates for 35-degree hyperboloid 




The pressure distribution along the wall is illustrated in figure 
52 for the three entry conditions. The results indicate that blunting 
of the body slightly increases the wall pressure for conditions of 
Z = 116 and 138 km. Results of the initial profile and the profile 
with uniform mass loss were so close that the difference could not 
be shown in figure 52. 

The results of radiative heating are presented in figure 53. 

These results follow the general trend of the sphere cone results. 

The differences in stagnation-point heating rates for the initial 
and blunted profiles are found to be 4.9, 4.1, and 3.2 percent for 
Z = 109, 116, and 138 km, respectively. Results of the profile with 
uniform mass loss are not seen to be significantly different from the 
results of the initial profile. 

Cqmparisoji of peak heating results . - Peak heating results for 
the 45-degree sphere cone, 35-degree hyperboloid, and 45-degree ellipsoi 
are compared in figures 54 to 58. Results for the 45-degree ellipsoid 
are seen to follow the same general trend as results for the 45-degree 
sphere cone; and, in comparison, results for the 35-degree hyperboloid 
are seen to exhibit a relatively smoother trend. In the stagnation 
region, all results for the ellipsoid are seen to be higher than 
results of the other two-body shapes. 

Variations in the shock-standoff distance, illustrated in figure 
54, indicate that the standoff distances for the ellipsoid are much 
greater than for the sphere cone, and, for the most part, the hyper- 
boloid results fall between these two results. For the ellipsoid, 
the effect of blunting is seen to be quite pronounced in the stagnation 
region. As such, one would expect a higher stagnation-region heating 
rate for the ellipsoid. 

The shock-temperature variations, illustrated in figure 55, 
indicate that the temperatures are higher for the ellipsoid near the 
stagnation region, but they fall between the results of the hyperboloid 
and sphere cone between s = 0.6 and 1.2. After location s = 1.2, 
the results are slightly lower than the results of the sphere cone. 
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Figure 55. Comparison of shock temperature for Z = 116 km. 








Figure 56. Comparison of shock density for Z = 116 
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In general, the shape change is seen to have greater effect on the 
temperature distribution for the ellipsoid than for the sphere cone. 
Variations in the shock density (fig. 56) follow essentially the same 
trend as the temperature variations. The effects of shape change 
for the ellipsoid, however, are not as pronounced as for the sphere 
cone. For the ellipsoid, blunting results in a maximum increase in 
density of 1.2 percent at s = 0.8. The results of uniform mass loss, 
however, do not show any significant change. 

The pressure distribution along the body surface is illustrated 
in figure 57. The results show that the shape change has a considerable 
effect on the pressure variations for the ellipsoid in the range from 
s = 0.6 to 1.2. The total drag for the ellipsoid, however, may not 
be greater than that for the sphere cone. As noted earlier, the shape 
change does not have significant influence on the pressure distri- 
bution for the hyperboloid, but the total drag for this shape can 
be higher than that for the other two shapes. 

The radiative heating rates for the three entry shapes are compared 
in figure 58. As expected, the radiative heating rates for the ellipsoi 
are comparatively higher in the stagnation region. In the downstream 
region, however, the results fall between the results of the hyperboloid 
and sphere cone. For the ellipsoid, blunting results in a maximum 
increase in heating of 7 percent at s = 0.8. In general, the increase 
in heating rates due to shape change is seen to be greater for the 
sphere cone and ellipsoid than for the hyperboloid. Also, the shape 
change is seen to have considerably more effect on heating of the 
afterbody for the sphere cone and ellipsoid than for the hyperboloid. 

The results further indicate that the total radiative heating load 
(i.e., the total radiative heat input) to the entry body will be 
comparatively higher for the ellipsoid, and this will be followed by 
the results for the hyperboloid and sphere cone, respectively. 

The results presented here indicate that uniform mass loss 
resulting in a shape that corresponds closely to the initial profile 
does not affect the shock-standoff distance, temperature, density, 




Figure 57. Comparison of surface pressure for Z = 116 km. 
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Figure 58. Comparison of radiative heating rates for Z = 116 km. 


and pressure distribution along the body appreciably. In such cases, 
radiative heating rates to the body are not different significantly 
from the heating rates to the body with initial configurations. 

Blunting of the nose region, however, is seen to have a significant 
influence on the entire flow phenomena in the stagnation and downstream 
regions. It is noted that the shock-standoff distance increases with 
increasing nose blunting. While nose blunting results in increasing 
the shock temperature all along the body, its influence on increasing 
the density is significant only in the stagnation region. In the 
downstream regions, the velocity, density, temperature, and pressure 
are altered significantly across the shock layer because of changes in 
the probe configurations. In most cases, considerable increase in 
radiative heating rates is noted in the stagnation as well as downstream 
regions due to severe nose blunting. Blunting of the entry body is 
seen to increase the wall pressure distribution significantly. But, 
its effect on stagnation-region and afterbody pressure distribution 
is relatively small. However, it is possible for the blunted sphere 
cone and ellipsoid to experience relatively higher total drag. 

Blunting is seen to increase the radiative heating rates all along the 
body for all configurations considered. But, in the stagnation region, 
the increase is relatively higher for the ellipsoid and sphere cone 
than for the hyperboloid. The shape change is seen to have considerably 
more effect on heating of the afterbody for the sphere cone and ellipsoid 
than for the hyperboloid. It is further noted that the total radiative 
heating load to the body will be comparatively higher for the ellipsoid 
followed by that for the hyperboloid and sphere cone. 
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CONCLUSIONS 


Results were obtained to investigate the influence of simplified 
radiation models, nonequilibrium radiative energy transfer, and probe 
condiguration changes on the flow properties and the heating rates in 
the stagnation and downstream regions of a Jovian entry body. Results 
obtained by using a simple, 30-step, radiation absorption model are 
in good agreement with results of other sophisticated models available 
in the literature. It is found that use of the present model reduces 
the computational time significantly. However, use of this simplified 
model is recommended only for general parametric studies. 

The radiative transfer equation has been formulated under the 
nonlocal thermodynamic equilibrium [NLTE) conditions. The NLTE 
effects are seen to enter through the absorption coefficient and the 
source function. The NLTE source function is expressed in terms of 
the Planck function, an NLTE parameter that measures the relative 
importance of the collisional and radiative deactivation processes 
in the gas and the influence factors arising from the higher level 
energy transitions. The Influence of NLTE on the entire shock-layer 
flow phenomena is investigated by neglecting the contributions of abla- 
tive products as well as by including them. 

The results obtained in the absence of the ablative species in 
the shock layer indicate that the NLTE effects are greater closer to 
the body than near the shock. The influence of NLTE, in general, 
is to reduce the convective and radiative heating to the entry body, 
and a significant reduction in radiative heating is noted. The NLTE 
effects are greater for higher entry conditions. The NLTE results, 
however, are not influenced by the change in the forebody configurations. 

The viscous shock- layer equations with coupled ablation and mass 
injection (for the entry probe with carbon phenolic heat shield) are 
solved under the NLTE conditions. The Swan (0,0), Freymark (0,1), 
and Mul liken (0,0) bands of the band systems are treated to be in 
nonequilibrium in the ablation layer. Flow-field results obtained 


for the peak-heating conditions (t = 111.3 sec) indicate that the 
temperature distribution in the shock layer is lower under NLTE conditions, 
tions. Similar behavior is also noticed for the enthalpy distribution. 

It is found that NLTE increases the density in the shock layer, but 
has no influence on the pressure variation. The radiative heating to 
the entry body is increased significantly because of NLTE; this, 
in turn, results in increased mass loss from the body. 

For investigating the influence of shape change of the entry probe 
on the flow field, different initial configurations (45-degree sphere 
cone, 35-degree hyperboloid, and 45-degree ellipsoid) for the entry 
probe were considered, and results were obtained for three different 
entry conditions (Z = 109. 116, and 138 km). The results indicate 
that uniform mass loss resulting in a shape that corresponds closely 
to the initial profile does not affect the shock-standoff distance, 
temperature, density, and pressure distribution along the body 
appreciably. Blunting of the nose, however, is seen to have a signi- 
ficant influence on the entire flow phenomena in the stagnation and 
downstream regions. Considerable increase in radiative heating rates 
is noted in the stagnation as well as downstream regions for all 
configurations considered. In the stagnation region, however, the 
increase is relatively higher for the ellipsoid and sphere cone than 
for the hyperboloid. It is concluded that, due. to the shape change, 
the total heating load to the body will be higher for the ellipsoid 
than the hyperboloid and sphere cone. 
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APPENDIX 


FINITE-DIFFERENCE SCHEME FOR VISCOUS 
RADIATING SHOCK- LAYER FLOW 


The solution of the second-order partial differential equation 
expressed by equation Cli3) is obtained by employing the implicit 
finite-difference scheme. For this purpose, the shock layer is 
considered as a network of nodal points with a variable grid space 
in the n-direction. The scheme is shown in figure A.l, where m 
is a station measured along the body surface and n denotes the 
station normal to the body surface. The derivatives are converted 
to finite-difference form by using Taylor’s series expansions. 

Thus, unequal space central difference equations in the n-direction 
at point m, n can be written as 
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A typical difference equation is obtained by substituting equations 
in equation (113) as 
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Figure A.l. 


Finite-difference representation of flow field 




where 


An = (2 aiAnn_p/[Arin(Ann * Ann.p] 

Bn = -[2 - aiCAn^ - ] /^AHn^nn.^ - a2 - (a4/A?j^_^)] 


C = (2 - aiAn )/[An i (An + An ,)] 
n ^ n-^ '■ n-1 ^ n n-1-^ 
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Now, if it is assumed that 


W = E W + F 

m,n n m,n+l n 

or 


(A. 3) 


W = E W + F , 

m,n-l n-1 m,n n-1 


(A. 4) 


then by substituting equation (A. 4) into equation (A. 2) there is obtained 


W = [-A /(B + C E J](W J 
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(A. 5) 


By comparing equations (A. 3) and (A. 5), one finds 


E = -A /(B - C E J 

n n n n n-1 


F = C-D - C F J/(B + C E ) 
n n n n-1 n n n-i 


(A. 6) 
(A. 7) 


Now, since E and F are known form the boundary conditions, E^ 

and F can be calculated from equations (A. 6) and (A. 7). The 
n 

quantities W at point m, n can now be calculated from 
m,n ^ 

equation (A. 3). 
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The overall solution procedure starts with evaluation of the flow 
properties iiranediately behind the shock by using the Rankine-Hugoniot 
relations. With known shock and body surface conditions, each of the 
second-order partial differential equations is integrated numerically 
by using the tridiagonal formalism of equation (113) and following the 
procedure described by equations (A. 2) to (A. 7). As mentioned before, 
the solutions are obtained first for the stagnation streamline. With 
this solution providing the initial conditions, the solution is marched 
downstream to the desired body location. The first solution pass 
provides only an approximate flow-field solution. This is because in 
the first solution pass the thin shock- layer form of the normal momentum 
equation is used, the stagnation streamline solution is assumed to be 
independent of downstream influence, dy^/d4 is equated to zero at each 
body station, and the shock angle a is assumed to be the same as the 
body angle 9. All these assumptions are removed by making additional 
solution passes. 

In the first solution pass, the viscous shock-layer equations are 
solved at any location m after obtaining the shock conditions from the 
free-stream conditions. The converged solutions at station (m-1) are 
used as the initial guess for the solution at station m. The solution 
is iterated locally until convergence is achieved. For the stagnation 
streamline, guess values for dependent variables are used to start the 
solution. In the first local iteration, (9n^/9^) and (9W/9C) are assumed 
to be zero. The energy equation is then integrated numerically to 
obtain a new temperature. By using this temperature, new values of 
thermodynamic and transport properties are calculated. Next, the x- 
momentum equation is integrated to find the u-component of velocity. The 
continuity equation is used to obtain both the shock standoff distance 
and the v-component of velocity. The pressure P is determined by 
integrating the normal momentum equation. Then the equation of state is 
used to determine the density value. 

With known stagnation streamline solution and body surface and 
shock conditions, the above procedure is used to find solutions for any 
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body location m. The flow chart for the computational procedure is 
shown in figures A. 2 and A. 3. Further details of this procedure and 
flow chart are given in reference 55. The flow chart for the NLTE 
radiation computation is shown in figure A. 4, and the integrals used 
in this chart are defined in figure A. 5. 
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Figure A. 2. Flow chart for shock- layer solution procedure 
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Figure A. 3. Flow chart for 






















Figure A. 4. Flow chart for subroutine NRAD. 
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Figure A. 4. (Concluded). 
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